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FOREWORD 


The Department of Physics and Astronomy of the University of 
Maryland and NASA Goddard Space Flight Center jointly organized a 
lecture series comprising an introductory course in high-energy 
astrophysics at the graduate level which was presented at the Uni- 
versity of Maryland in the summer of 1968. Throughout the lectures, 
notes were distributed to the students to serve as basic text material 
for the course. This publication is an expanded and edited version of 
those notes. 
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I. CURRENT PROBLEMS OF HIGH-ENERGY ASTROPHYSICS 


H. Ogelman 

NASA Goddard Space Flight Center 
Greenbelt, Maryland 


INTRODUCTION 

If we drop the prefix "astro," high-energy physics may be de- 
fined as that branch of physics concerned with the interactions be- 
tween particles having relative energies in excess of a few hundred 
MeV, the energy threshold at which the familiar particles of nuclear 
physics, such as the proton and neutron, start yielding pions, kaons, 
etc. Since the fundamental interactions that produce these elementary 
particles seem to be distinct from those observed at lower energies, 
and since the well-controlled laboratory environment makes it pos- 
sible to focus attention on a single process, such a categorization is 
justified; it aids the scientists who can best understand a phenomenon 
by restricting the scope of their interests. 

In astrophysics, however, where the goal is to describe the 
physics of the universe, the crucial advantage of controlled experi- 
mentation is not available. In order to increase our understanding 
of astrophysics, we must employ information from all available 
channels; it would be harmful to restrict ourselves to a particular 
energy band or fundamental process. In defense of specialization, 
however, it is reasonable to assume that certain physical processes 
will achieve prominence in certain regions of the universe. Concen- 
trating one's attention on a particular type of process and then re- 
lating one's findings to those of colleagues specializing in other types 
of process can therefore be expected to provide fruitful results. 

The field of high-energy astrophysics, consequently, covers proc- 
esses involving energies of 10 3 eV to 10 20 eV — 17 decades of energy — 
as compared to the 4 decades of energy covered by present-day high- 
energy physics (10 7 to 10 11 eV). 

Historically, high-energy astrophysics has developed from cosmic 
ray physics. Cosmic rays were discovered over 60 years ago by re- 
searchers studying the discharge of carefully insulated electroscopes. 
When Hess flew an electroscope in a balloon in 1911 and showed that 
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the discharge rate increased as a function of altitude, it became evi- 
dent that the cause of the discharge was extraterrestrial. During the 
past 50 years, extensive work in the field has revealed greatly the 
astrophysical significance of cosmic rays. The secondaries of cosmic 
rays were actually the starting point of elementary particle physics; 
particles such as muons, positrons, pions, kaons, and some hyperons 
were first discovered in cosmic rays. 

With the advances in high-energy accelerators during the past 
few decades, the aspirations of cosmic-ray studies have shifted from 
elementary particle physics to astrophysics. However, above 10 1 1 eV, 
cosmic rays are still the only accessible high-energy beam; they are, 
therefore, still contributing to the field of high-energy physics, in 
particular to the search for theoretically predicted particles such as 
quarks, intermediate bosons, and monopoles. 

Cosmic-ray energies currently under investigation in high-energy 
astrophysics range from 10 s to 10 2 ° eV. The universe, however, is 
a complicated laboratory system in which the beam, the target, and 
the secondaries produced in the reactions do not always exhibit the 
patterns with which we are familiar in the laboratory. For example, 
a seemingly insignificant 10" 3 eV microwave photon is as potent to 
a 10 19 eV cosmic ray proton as a 300 MeV photop produced by an ac- 
celerator is to a proton in a hydrogen bubble chamber. In other 
cases, the cosmic target may be the weak magnetic fields of the 
galaxy, where an incident electron of 10 10 eV could, by interacting 
with the field, produce a radio photon of 10" 7 eV. The information 
relevant to the high-energy phenomena could cover the energy inter- 
val from 10 7 to 10 2 ° eV~26 decades of energy! 

What types of channels would be able to link us to the rest of the 
universe and to convey valuable information? Among the discovered 
elementary particles, we must look for stable ones that will not de- 
cay in flight. If a particle has a natural lifetime r, it can carry in- 
formation up to a distance R, given by 

R <%/ 

- ycT , 


where y is the Lorentz factor of the particle [l -• ( v/c) 2 ]" 1 72 , v is 
the particle velocity, and c is the velocity of light. For the extreme 
case of 10 20 eV protons, which is on the verge of present-day detection 
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schemes, y is 10 11 . Assuming that this is a practical limitation on 
currently observable values, and that we are interested in reaching 
the limits of our galaxy (i.e., R -10 22 cm), this gives 


t > 3 sec . 


Among the measured elementary particle lifetimes, the value falls 
between 10 3 sec for the neutron and 2 x 10" 6 sec for the muon, justi- 
fying the previous crude argument. Effectively then, we can eliminate 
the mesons, hyperons, and muons as a part of the incident primary 
radiation because of their short lifetimes. The rest of the possible 
primaries are the photon, electron, positron, neutrinos, proton, and, 
to a certain extent, the neutron. Table 1 summarizes various rele- 
vant parameters of these particles. 


Table 1. Summary of Relevant Parameters 
of Photons, Leptons, and Nucleons 


Particle 

Symbol 

Spin 

Rest mass 
(MeV) 

Mean life 
(sec) 

Decay mode 

Photons 

7 

1 

0 

CO 



Leptons 

e~ , e + 

1/2 

0.51 

CO 

— 



1/2 

0 

CO 

— 



1/2 

0 

CO 

— 

Nucleons 

p 

1/2 

938.2 

CO 

— 


n 

1/2 

939.5 

1000 

n -* e~ + v + p 
e ^ 


Besides these elementary particles, any stable or quasi-stable 
(in the sense of the neutron) combination of elementary particles, 
such as alpha particles or heavy nuclei, should be observable in the 
primary beam of cosmic rays. Further coalescence of particles to 
atoms, atoms to molecules, molecules to grains, and so on, is pos- 
sible and, if spared by destructive collisions, these objects may also 
be present among the primaries. 
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Another interesting possible source of information could be the 
graviton, the stable quantum of the gravitational field. Although no 
graviton has been observed in a quantum state, its manifestations in 
the form of gravitational force field are well known. Pioneering 
astrophysical research in the field of gravitational waves is presently 
in process (References 1 and 2). 

A review of existing information on primary cosmic radiation and 
a discussion of the general features of our galaxy and the metagalaxy 
follow. 

PRIMARY COSMIC RAYS 

The purpose of measuring the primary cosmic-ray flux is to dis- 
cover the intensity I of each component (proton, electron, photon, etc.) 
of the flux as a function of energy E, direction 0,<P, time t, distance R, 
and spina, (commonly known as polarization). Unfortunately, tech- 
nological limitations do not allow measurement of all of these quanti- 
ties. The difficulties arise mainly for the following reasons: 

1. The earth's atmosphere is many interaction lengths for the 
nucleonic component (13 mean free pathlengths for protons) and the 
high-energy electromagnetic component (25 mean free pathlengths at 
high energies). This fact forces experimenters to construct apparatus 
that can be exposed at the top of the atmosphere via balloon, rocket, 
or satellite, thus limiting the size. These inconveniences currently 
limit the possible flux measurements to values exceeding 10" 8 par- 
ticles per cm 2 -sec-sr. 

2. The very low flux of certain components relative to the proton 
component makes it difficult to distinguish them from the secondaries 
produced in the atmosphere or the apparatus. 

3. The solar wind, with its complicated magnetic fields, affects 
those charged primaries having rigidities below a few GeV and pre- 
vents observation of the true galactic distribution. 

Despite these difficulties, a wealth of data about primary cosmic 
rays has been acquired. Some of these data are summarized below. 
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Protons and Nuclei 

The most abundant component of cosmic rays is the proton. Be- 
sides the protons, cosmic rays contain 6 percent helium nuclei and 
1 percent higher Z nucleons. The direct composition measurements 
cover the energy range of 10 8 eV per nucleon to about 10 10 eV per 
nucleon. Table 1 of Chapter II summarizes the relative intensities of 
the various cosmic-ray components and gives cosmic abundance esti- 
mates, The excess of heavy elements in cosmic rays relative to 
cosmic abundances is generally interpreted as being caused by ac- 
celerating regions, and the excess of light elements as being caused 
by fragmentation of heavier elements during their passage through 
ambient gas. 

The spectral form of the incident cosmic radiation can be 
expressed as a power law (Reference 3). In the energy range 
10 10 eV < E <10 15 eV, it takes the form 


I(>E) 


KE~ ( 1 ) 



where E is measured in Ev. In the range 10 1S eV < E < 10 16 eV, 


I( >E) 



and, for E > 10 16 eV, 


I(>E) 



the total number of particles being 1500/m 2 -sec-sr. The primary 
spectrum is shown in Figure 1. The breaks occurring in the spectrum 
at 10 15 and 10 18 eV are generally interpreted as being caused by the 
escape of the galactic component and the penetration of the extra- 
galactic component, respectively (Reference 4). 



LOG E (eV) 


Figure 1- Integral primary energy spectrum (taken from 
Figure 1 of Reference 4). 
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Because of the sharp drop in intensity as a function of energy, 
direct measurements have been carried out only to 10 14 eV. Above 
this value, extensive air showers have provided most of the informa- 
tion (Reference 5). At these energies, the atmosphere becomes a 
part of the detector system and makes it feasible to extend the energy 
measurements out to 10 20 eV. Figure 2 shows schematically the de- 
velopment of an extensive air shower in the atmosphere. Typically, 
the incident nucleon interacts with air molecules and produces pions 
which in turn decay into electromagnetic channels via the tt ± -* /x* - e* 
or the 77° - 2y processes. Progressively, as the average energy per 
particle decreases, their numbers increase, reach a maximum, and 
then die out. During this development, the particles form a coherent 
front recognizable by detectors placed on the ground. 

Two advantages of the atmosphere as a detector are evident. One 
is its extensive depth, which allows even 10 20 eV showers to develop 
well within the atmosphere; the other is its low density, which allows 
the spreading of the shower particles, enabling a small number of 
detectors to cover effectively a large area. A primary proton of 10 19 
eV produces an easily detectable flux of secondaries 1 km from its 
impact point. Since a single detector would not be able to distinguish 
the difference between a large shower far away and a smaller one 


PRIMARY 

PARTICLE 



PHOTONS 


Figure 2— Schematic diagram of the propagation of an 
extensive air shower (Wolfendale, 1963). 
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close by, the technique used is to distribute a number of detectors 
over an area of 10 or more square kilometers and look for coincident 
pulses. Moreover, the relative timing of the pulses in the shower 
front allows the determination of the primary direction. The density 
of particles at various detectors allows estimation of the total num- 
ber of secondaries N which can then be related to the energy of the 
primary particle by theoretical arguments. 

The penalty paid for the convenience of air-shower detection is 
the loss of some information regarding the nature of the primary. 
Still, some signature of the primary is retained. For example, 
showers with relatively small numbers of muons could be attributed 
to an incident electron or photon (References 7 and 8) since, in such 
a shower, electromagnetic cascade development would be enhanced 
relative to nuclear cascade development. A shower initiated by a 
multi-nucleon particle might show separated peaks of electron den- 
sity near its core because of individual showers developed by each 
nucleon (Reference 9). 

The results of air- shower measurements have indicated the 
plausibility of a significant flux of protons and high-Z primaries of 
around 10 14 to 10 15 eV total energies; however, the interpretation of 
the data is complex, and the results are not well confirmed. 

The isotropy of cosmic rays in the energy ranges discussed has 
been remarkable. The usual definition of isotropy 8 is 


„ ^raax ^min 

8 “ I + I . 

max min 


The conservative upper limits of 8 at various energies E are 
(Reference 10): 

8 < 0. 1% , E < 10 14 eV ; 

8 < 1% , E ~ 10 16 eV ; 

S < 3% , E ~ 10 17 eV . 
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However, there are some indications that around 10 12 eV there 
is a small anisotropy, of fractions of a percent in amplitude, from the 
direction of the galactic center. When one considers that, in the weak 
magnetic fields of about 3 x 10 -6 gauss present in our galaxy, the 
radius of curvature of a primary with energy 10 19 eV is the same 
order of magnitude as the size of our galaxy, the difficulty of accelera- 
tion and confining these particles within our galaxy becomes under- 
standable. 

Another topic under investigation is the long-time scale variation 
of cosmic radiation. The basic measurements have been conducted 
on meteorites which have been traveling in space for long periods of 
time. If one measures the induced radioactivity for different elements 
having various decay lifetimes, it is possible to interpret the concen- 
tration as induced by cosmic-ray bombardment. Measurements of 
decay half-lives from 12.5 years to 10 9 years indicate that cosmic 
rays have been bombarding the solar system at about their present 
rate for over a billion years. 

For a more complete discussion of this whole area, see Chapter 
II, "Origin, Composition, and Propagation of the Nuclear Component 
of Cosmic Rays,” by C. Fichtel. 


Electrons 

Significant measurements have been made during the past several 
years, resulting in the determination of the primary electron spectrum 
in the energy range 10 6 to 10 11 eV (References 11, 12, and 13). It is 
generally agreed that the flux below a few GeV is dominated by solar 
effects. The higher energy flux, believed to be galactic in origin, can 
be represented by a power law: 

( £ \- 2 . 3 ± 0.2 

^ j electrons m “ 2 - sec" 1 - sr" 1 -GeV " 1 . 


Figure 3 shows the measured spectrum on top of the atmosphere. 
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Figure 3-The differential energy spectrum of primary electrons. Re- 
sults of experiments are indicated; Minnesota (Webber and Chotkowski 
1966), Chicago (L'Heureux 1968); Leiden (Blecker et al. 1967); and 
Bombay (Daniel and Stephens 1966). Taken from Scheepmaher and 
Tanaka, 1968. 
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The electron-to-positron ratios have also been measured at 
various energies; Table 2 (from Reference 14) summarizes these 
measurements. 


Table 2. Summary of Measurements of 
Electron-to-Positron Ratios 


Electron energy 
(eV) 

e + /(e + + e ) 

1.75 x 10 8 
6 x lo 8 
2 x 10 9 

0.43 ± 0.21 
0.30 ± 0.14 
0.33 ± 0.16 


Independent measurements on the electron spectrum can be made 
from radio astronomy data by interpreting the galactic background 
noise as being caused by electron synchrotron radiation. Presently, 
a consistent picture is emerging, yielding information on the lifetimes 
of these electrons in the galaxy and on their production mechanisms. 
For a more detailed treatment of this subject, see Chapter III: 
"Cosmic-Ray Electrons and Related Astrophysical Problems," by 
R. Ramaty. 

Photons 

In contrast to the charged particles, which are deflected by in- 
terstellar magnetic fields , photons travel in straight lines and can 
thus be traced to their origin. The visible photons are very aniso- 
tropic, exhibiting their points of origin (stars, nebulae, galaxies, 
etc.). The radio-frequency photons exhibit similar anisotrophy. 

Since the atmosphere is transparent to electromagnetic radiation 
between 10' 8 and 10“ 3 eV and around 1 eV, ground-based astronomy, 
both optical and radio, has yielded much information. There are 
several outstanding features of this low-energy photon information 
that are related to high-energy astrophysics. The measurement of 
the anisotropic radio waves and their interpretation as synchrotron 
radiation emitted by high-energy electrons in the weak magnetic 
fields of the sources has been one of the most significant achievements 
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of high-energy astrophysics. The rate of energy emission by an 
electron of energy E in a magnetic field His given by 


dE 

dt 


- 3.8 x 10" 6 H 2 ( gauss )E 2 (GeV) 10 9 eV/sec , 


and the frequency of emission is strongly peaked around: 


v c - 1.6 x 10 13 H( gauss) E 2 (GeV)Hz . 


Historically, it was the Crab Nebula, a supernova remnant over 
900 years old, to which these measurements and theoretical calcula- 
tions were first applied. Indirect arguments concerning the energy 
balance of the gas and the magnetic pressure lead to magnetic fields 
of 10” 3 to 10” 4 gauss. These values of H imply that 10 9 eV electrons 
emit the radio noise and that 10 u to 10 12 eV electrons emit the blue 
continium light received from the Crab Nebula. The synchrotron 
emission mechanism is identified primarily by the polarization of 
the radiation. Experimental measurements of this polarization have 
confirmed the synchrotron radiation theory. 

The same type of measurements, applied to other discrete radio- 
frequency emitting regions exhibiting nonthermal spectra, proved the 
existence of high-energy electrons in various other localities in our 
galaxy and in other galaxies. Up to the present time, radio measure- 
ments have been the only means of localizing the sources of high- 
energy particles. However, radio emission, being a secondary effect, 
is connected with important parameters such as the energy distribu- 
tion of the electrons through the uncertain estimates of local magnetic 
field strengths. It is, of first rate importance to observe the high- 
energy flux itself and to associate it with the sources. Recent meas- 
urements of gamma rays above 100 MeV by the OSO-in satellite-borne 
detector indicate a strong anisotropy towards the galactic center, 
thereby constituting the first clear detection of high-energy photons 
(Reference 15). 

Another recent discovery having an important bearing on high- 
energy particle interactions is the existence of 3°K black body 
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radiation, where the peak of the radiation is at about 7 x 10" 4 eV 
(Reference 16). This black body radiation is generally interpreted 
as being universal and a consequence of the decoupling of matter and 
radiation at an early epoch of the universe (Reference 17). Because 
of its high photon density (600/cm 3 ), this radiation renders the uni- 
verse opaque to particles with which it has resonant interactions 
(Reference 18). 

The diffuse photons of lower energy pervading the universe are 
shown in Figure 4. In our galaxy, the optical and lower frequency 
radio noise is enhanced relative to the microwave noise. However, 
the photon spectrum is still divided into four distinct energy regions, 
namely nonthermal radio, microwave, optical, and X-rays. For a 
detailed treatment of this topic see Chapter IV, "Cosmic-Ray Photons, 
by E. Boldt. 



Figure 4-Power spectrum of metagalactic photons as assessed 
by Gould and Schreder, 1967. 
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Neutrinos 

Being a neutral particle, the neutrino should offer the same di- 
rectional advantages as the photon. However, it exhibits weak inter- 
actions only. The total charged-particle production cross-sections 
are of the order 10“ 3 8 E„ cm 2 per neutrino, where E v is the neutrino 
energy in GeV (Reference 5). This is the theoretical cross-section 
inferred from Fermi's universal weak interaction theory and is valid 
for energies in excess of 1 GeV; at lower energies, the cross-section 
is proportional to E v 2 . The magnitude of the problem is obvious; kiloton 
detectors would be necessary to detect such particles. The estimates 
of the neutrino flux.from interstellar space at energies above 1 GeV 
are similar to those of gamma rays, amounting to 10’ 4 neutrinos 
cm -2 -sec _1 -sr -1 . However, the atmospheric secondaries should be 
as intense as 10" 1 neutrino cm" 2 -sec" 1 -sr" 1 . In view of this dif- 
ficulty, high-energy neutrino astronomy could only bear fruit if 
kiloton detectors could be assembled in space or on the moon. High- 
energy neutrinos have been observed over a mile underground, but, 
these are generally considered to be atmospheric secondaries (Ref- 
erence 20). 

Low-energy solar neutrinos that are produced by the ordinary 
nuclear reactions involved in hydrogen burning have been a topic of 
recent interest (Reference 21). Detectors have been built that could 
detect the Ar 37 formed by the reaction 


Cl 37 + v - Ar 37 + e" 


The negative results of these experiments challenge the theoretical 
predictions (Reference 22). 


THE ASTRONOMICAL SETTING 

Our sun, with the parameters 


= 

2 x 

10 33 g, 

^ = 

4 x 

10 33 ergs/sec 



CURRENT PROBLEMS 


15 


and 


R 0 - 7 x 10 10 cm , 

is one of the 10 11 stars filling our galaxy. A schematic drawing of the 
galaxy, indicating its dimensions and shape, is shown in Figure 5. The 
star distribution in the galaxy is concentrated in a disk about 300 pc 
thick (1 pc = 3.1 x 10 18 cm) with an average density of 0.1 star/pc 3 . 
The density of the stars in the nucleus seems to increase to 100 stars/ 
pc 3 . Most of the stars in the disk and the nucleus belong to population 
I, which is the younger type of star. A small fraction of the stars 
form globular clusters which are located away from the disk. These 
stars are generally of population II, i.e. , older stars. 



Figure 5-Schematic representation of our galaxy 
(clotted lines show the extent of proposed halo). 
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The spiral-arm structure of our galaxy is more evident when 
one examines the interstellar gas distribution. Studies of the gas 
distribution became possible after the discovery of the monochromatic 
21-cm hyperfine transition line emitted by the neutral hydrogen atom. 
Figure 6 shows the profiles of the 21-cm hydrogen line. Interpreting 
the shift in frequency as caused by galactic rotation, one can determine 
from the shape of the observed emission the neutral hydrogen density 
at that direction and distance. The resulting picture is shown in Fig- 
ure 7, clearly indicating the spiral structure of our galaxy. These 
regions of neutral hydrogen, called H t clouds, have average densities 
of 10 atoms/cm 3 and are at a temperature of 100°K. The random 
cloud velocities are of the order of 10 km/sec. The clouds seem to 
be lumpy with scales of 10 pc. The neutral hydrogen occupies about 
10 percent of the volume of the disc, the other 90 percent being 
ionized hydrogen clouds called H n regions. In these H n regions, 
the density is estimated to be ^0.2 ion/cm 3 , and the temperature to 
be ~10 4 °K. The existence of the molecular hydrogen is still an open 
question. Published estimates vary from 10 percent of the atomic 
hydrogen concentration to 10 times this concentration. If low-energy 
cosmic rays ('TOO MeV) pervade the galaxy, it is extremely difficult 
to understand how the molecules could be preserved against the dis- 
ruptive ionizing effects of the radiation (Reference 25). 

The other constituent of the interstellar medium is the grains or 
macromolecules which cause the observed interstellar reddening. 

The average density of these is estimated to be less than 10' 2 5 gm/cm 3 . 

The existence of a spherical halo surrounding our galaxy was 
first suspected as a consequence of the slow decrease of radio in- 
tensity in directions away from the disk. Although some of the other 
spiral galaxies exhibit such a halo, others do not. At the present, it 
is regarded as a theoretically attractive characteristic of our galaxy 
although experimental evidence is inconclusive. The plausible shape 
and size of the halo is a spheroid with a major radius of 15 kpc and 
a minor radius of 8 kpc. The density of hydrogen outside the disk is 
estimated to be 10' 2 6 gm/cm 3 or less. 

The nucleus of the galaxy has a diameter of approximately 8 pc 
as measured by thermal radio emission of ionized hydrogen, the 
concentration of which may be as high as 10 3 atoms per cubic 
centimeter. 
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Figure 7— Contour map of neutral-hydrogen density in the galactic system 

(Reference 24). 


Magnetic fields in the galaxy were first postulated by Fermi in 
1949 to explain the isotropy of cosmic rays (Reference 26). Today we 
have independent evidence in the polarization of distant starlight, in- 
terpreted as the preferential absorption of grains oriented by the 
magnetic fields, and nonthermal radio emission, interpreted as syn- 
chrotron radiation of high-energy electrons in galactic magnetic 
fields. Measurements of these quantities and theoretical energy 
balance arguments seem to favor an average field of 3 x 10' 6 gauss, 
with field lines running in the plane of the galaxy. It seems plausible 
that, in the spiral arms, the field strength may go up to 10' 5 gauss. 
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The photon densities in our galaxy can be listed in the radio, 
microwave, and optical regions. In the disk, optical photons have a 
density of 0.2 photon/cm 3 ; radio photons, about 0.1 photon/cm 3 . In 
the halo, these densities may be less by a factor of 3. Microwave 
radiation, being universal, pervades the entire galaxy at a density of 
600 photons/cm 3 * . 

In the vicinity of our galaxy, we find 15 other galaxies, forming 
with ours what is termed the local group, including the Magellanic 
Clouds and the well known Andromeda nebula. The size of this sys- 
tem is approximately 1 megaparsec. The average separation of these 
galaxies is approximately 50 times their diameter. 

Other clusters of galaxies, such as Virgo, have been observed to 
contain as many as 1000 galaxies. The other observed galaxies seem 
to recede from us at a rate V proportional to their distance R: 


V = HR , 

where H is the Hubble constant, 10 8 cm/sec-Mpc. Extrapolating this 
relationship to a velocity equal to that of light gives 10 2 8 cm as the 
observable limit to the universe. The galaxies observed within this 
volume show a variety of different properties and power outputs. 

For example, it is believed that quasars are distant, relatively small 
galaxies with energy outputs of 10 5 2 ergs/sec as compared with an 
average galaxy's output of 10 45 ergs/sec. It is generally believed 
that, in these peculiar galaxies, high-energy processes much more 
intense than those encountered in our own galaxy are dominant. 
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II. ORIGIN, COMPOSITION, AND PROPAGATION OF THE 
NUCLEAR COMPONENT OF COSMIC RAYS 

C. Fichte 1 

NASA Goddard Space Flight Center 
Greenbelt, Maryland 

NUCLEAR COMPOSITION 

In studying the origin and propagation of the cosmic radiation, it 
is valuable to begin by reviewing the composition of the nucleonic com- 
ponent. Table 1 gives the relative abundances of some of the nuclear 
components of interest here, as well as some relevant ratios. 


Table 1. Relative Abundances of Indicated Nuclear Components 
With a Base of 100 for He; and Ratios of Some Nuclear Groups 


Nucleus 

- 

High energy 
cosmic rays 

Universal abundances 

Protons 

(1.5 ± 0.2) x 10 3 

(0.7 to 1.3) x lo 3 

He 

100 

100 ± 10 

Li, Be, B 

1.6 ± 0.3 

10" 5 

C, N, 0, F 

6.7 ± 1.6 

1.0 

H(Z > 10) 

2.5 ± 0.3 

(0.1 to 0.4) 

VH(Z > 20) 

0.7 ± 0.2 

(0.005 to 0.02) 

He 3 /(He 3 + He 4 ) 

-0.2 

- 

Carbon/ Oxygen 

1.1 ± 0.1 

(0.1 to 0.4) 

(Z > 30)/(Z > 20) 

-0.4 x lo -3 

(0.4 to 1.3) x 10" 3 

(Z >40)/(Z > 20) 

-0.3 x IQ -4 

(0.5 to 3) x 10“ 4 


The Li, Be, B, and He 3 in the cosmic radiation are thought to be 
secondaries formed in interactions between cosmic ray nuclei and 
interstellar material (presumably mostly hydrogen). Although there 
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are some difficulties in explaining the amount of these nuclei as a 
function of energy, which will be discussed in the last part of this 
lecture, the abundances of these secondary nuclei can be explained in 
general terms by assuming that on the average the cosmic radiation 
has gone through about 4 ± 1 g/cm 2 of material. This somewhat un- 
usual unit, i.e. "g/cm 2 ," arises because it is the product of a length 
times a density, a quantity that determines the degree to which the 
heavier cosmic ray nuclei break up in interactions. With this figure 
and the current best estimate of the density of interstellar material in 
the galactic disk— about 1 atom/cm 3 — the average lifetime of the cos- 
mic radiation in the disk is found to be a few million years. This num- 
ber, together with our knowledge of the cosmic ray energy spectrum 
and flux, will tell us the rate at which cosmic rays must be supplied. 

If the cosmic radiation pervades the whole galaxy and not just the disk, 
it still has to be supplied at about the same rate because its lifetime 
is larger by approximately the same factor as the volume. The subject 
of the region of containment of the cosmic rays will be discussed later 
in this chapter. 

The volume of the disk will be taken as 

77 (l . 5 x 10 4 x 3 x 10 18 ) 2 x (400 x 3xl0 18 ) = 8 x 10 66 cm 3 . 


If there are 2 x. lo -10 cosmic rays per cubic centimeter, the total 
number of cosmic rays in the disk is a little more than 10 s7 and their 
energy is 8 x 10 66 eV or about 10 5 5 ergs. Using the lifetime derived 
above, the rate at which cosmic rays must be supplied is about 2 x 10 43 / 
sec, and the rate at which energy must be supplied is about 10 41 ergs/sec. 

The numbers quoted in the last paragraph refer to cosmic rays 
with energies above about 0.3 GeV/nucleon. At lower energies the 
cosmic radiation is strongly modulated by the outward flowing hot 
solar plasma called the "solar Wind." This is a problem in itself and 
will not be pursued here. The degree of modulation is the subject of 
considerable controversy, and estimates of the number of low energy 
cosmic rays outside the solar system vary by more than an order of 
magnitude. Fortunately, because their energy is so much lower than 
the average cosmic ray energy, the total cosmic ray energy is probably 
not increased by more than a factor of two even under the most ex- 
treme assumptions, and probably less. 
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Before leaving the subject of secondaries in the cosmic radiation, 
which is what led to the discussion just completed, a few other features 
will be mentioned, A significant fraction of the nuclei with charges 
just below iron are probably also secondaries on the basis of their un- 
usual abundances. Many of the odd nuclei from Z = 7 to 19 would also 
be expected to be secondaries, and the carbon-to-oxygen ratio should 
be enhanced, as indeed it appears to be. 

The excess of heavy nuclei in the cosmic radiation, relative to 
universal abundances, is seen to increase generally from the medium 
nuclei (where the excess relative to helium is about a factor of 6) to 
the iron group (where it is about a factor of 100). This increase is 
usually attributed to the character of the source region. The nuclei 
with charges greater than iron also show an excess relative to helium 
in the cosmic radiation, but the ratios of nuclei with Z > 30 and Z > 40 
relative to the iron group appear to be the same as the universal 
abundances. This feature is also thought to be a characteristic of the 
source region, and theoretical attempts to explain this composition 
are a subject of considerable current interest. 

Another feature of the nuclear component which enters the picture 
is that, except for small variations, the relative abundances of the 
various nuclei seem to vary little with energy from energies as low 
as 0.05 GeV/nucleon up to at least 50 GeV/nucleon, and there is some 
evidence that the relative abundances of the major groups do not 
change substantially even up to energies as high as 10 6 GeV/nucleon. 


ORIGIN 

Region of Containment 

There are basically three major regions in which the cosmic rays 
could be contained and be presumed to originate: a small region sur- 
rounding the solar system; bur own galaxy; and the universe. 

The solar system can be rejected for numerous reasons, including 
the following: 

1. The magnetic fields of the solar system cannot contain the 

highest energy particles, that is those above about 10 1 2 
GeV/nucleon. 
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2. The cosmic rays are most intense when the sun is least active. 

3. The composition of the cosmic radiation differs greatly from 
that of the sun and that of energetic solar particles (average 
energy between a few and tens of MeV/nucleon) which the sun 
is known to emit in association with some big solar flares. 

4. The energy of the particles in the solar particle events is 
very small compared to the average cosmic ray energy, and 
there is no known way that the sun could create such high 
energy particles and retain its present characteristics. Even 
if the sun could emit the higher energy particles, they would 
be highly anisotropic, since they would be virtually undeflected 
by the magnetic fields in the solar system. 

The universe as a whole seems unlikely as the region in which 
cosmic rays are localized for several reasons. 

Gamma-ray studies indicate that the density in interstellar space 
would have to be less than 10” 7 atom/cm 3 if the cosmic ray intensity 
were the same as in the vicinity of the solar system. This density 
seems too small to be consistent with x-ray measurements and other 
considerations but the question is not finally resolved. If cosmic rays 
are extragalactic in origin, the galactic magnetic field must connect 
to the metagalactic region to let cosmic rays into the galaxy and out 
again in about 10 6 years; but the nature of this connection is difficult 
to visualize since the galaxy has rotated about twenty-five times in 
the life of the universe. There is also the question of finding a source 
or sources to supply the needed energy. 

These and similar difficulties are strong objections, but they do 
not absolutely eliminate the universe or some local region of it— most 
probably the region in which the highest energy cosmic rays (E ~ 3 
x 10 18 eV/nucleon) are contained. For the present discussion, how- 
ever, the universe will also be rejected as the region of containment 
for the majority of cosmic rays. 

Turning again to the galaxy, two basic conditions must be satisfied. 
First, the cosmic rays must be produced at the desired rate with the 
necessary characteristics; and secondly they must be confined satis- 
factorily to the galaxy. With regard to the second condition, we may 
assume that the maximum rigidity that a particle can have and still 
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be contained easily in the galaxy is one which is an order of magni- 
tude smaller than that corresponding to a radius of gyration equal to 
the radius of the galaxy. This rigidity is 

R = { 300 r gal. ( Cm > H gal. (gauss)} • {^} 

Using r gal = 1.5 x 10 4 parsecs and H gal = 3 x 10" 6 gauss (this 
number is very poorly known), gives a rigidity of 5 x lo 18 eV; thus, 
all but particles with the very highest rigidity could be contained in 
the galaxy. Similar arguments would indicate that a particle with 
rigidity equal to about 3 x lo 16 eV or less could be contained in the 
disk. 

In examining the problem of containing cosmic rays in the disk, 
the various contributing pressures for expansion must be considered 
and their resultant gradient balanced against pg, the local density 
times the local value for the acceleration of gravity. The expansive 
pressures are those due to the cosmic ray gas, to the magnetic field, 
and to the kinetic motion of the interstellar particles . Parker has 
shown that for a magnetic field of (3 to 5) x lo" 6 gauss (estimated 
from Faraday rotation measurements and polarization measurements, 
assuming graphite grains are primarily responsible for the polariza- 
tion) , a balance occurs for an interstellar density of 1 to 2 hydrogen 
atoms/cm 3 . This is in good agreement with current estimates. If a 
similar calculation is made for the galactic halo, the relatively low 
density of material thought to exist there implies that the cosmic ray 
density is quite low; so, in effect, cosmic rays are not stored there if 
current ideas are correct. 

The magnetic field configuration in the disk is actually unstable 
and forms humps from which the cosmic rays can escape more easily. 
The cosmic rays must escape somehow, because they are being sup- 
plied continually and yet their average energy density remains con- 
stant. The unstable motion of the magnetic fields in the disk implies 
a filamentary structure for the galactic field, and there is some 
evidence to support that hypothesis. It is thought that the diffusion of 
cosmic rays in this field is rather slow, and that as a result the cosmic 
rays are a relatively local phenomenon. It is possible that the cosmic 
rays seen near the earth come mostly from sources in the galactic 
spiral arm segment in which our sun is located. The cosmic ray 
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bulk motion perpendicular to the disk is very slow in this picture, 
since the main component of the field is parallel to the plane of the 
disk, and the net motion is due as much to the pressure balance men- 
tioned above as to the result of scattering and spiraling in the magnetic 
field. 

Source 

Assuming then that the cosmic rays can be contained in the galaxy, 
what is their source? Current theories seem to favor supernovae be- 
cause they seem to be the only objects capable of supplying the neces- 
sary energy. From the theory of supernovae and from the radio ob- 
servations interpreted in terms of synchrotron emission of relativistic 
electrons, it seems reasonable that each supernova will supply from 
10 49 to 10 52 ergs to cosmic rays. Assuming a supernovae occurs 
once every 50 years in our galaxy, supernovae would be able to sup- 
ply cosmic ray energy at the rate of about 10 40 to 10 4 3 ergs/sec on 
the average. The energy supply rate calculated earlier in this chapter, 
10 41 ergs/sec, is within this range. Supernovae are also believed to 
be rich in heavy elements, and this is consistent with cosmic ray com- 
position. In order to proceed further, and in particular to consider 
both the problems of energy spectra and of the maximum energy, it is 
necessary to examine specific models for supernovae. 

Present theories of supernovae are based on the concept that a 
supernova explosion is initiated by a dynamic implosion of a massive 
star. This implosion results from the gravitational instability of the 
star, at the end of nuclear synthesis, wherein there is thermal de- 
composition of iron back into helium. Although there is some debate, 
the most popular current view is that the gravitational instability oc- 
curs because the neutrinos emitted in the electron capture process 
leading to neutron- rich material remove energy faster than quasi- 
static contraction can supply it. As a result, a dramatic implosion 
occurs and continues in virtual free-fall until the pressure of the 
Fermi gas neutrons in the core becomes high enough to stop it. 

A large gravitational energy is thus freed and transferred to the 
neutrinos. The neutrinos, whose mean free paths are smaller than 
the radius of the star, deposit their energy inside the star, and a 
shock wave is formed. The velocity of this shock wave increases 
towards the surface and becomes relativistic for a small fraction of 
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the outer envelope, which is then converted into cosmic rays. To 
discuss this process in detail is far beyond the scope of this lecture, 
but the results predict a power law spectrum in energy similar to 
that actually observed, and it is supposedly possible to obtain the 
highest energies by this process. The excess of heavy nuclei can 
also be explained by this model, but there is some question as to 
whether it predicts the correct composition in detail. 

One consequence of the hydromagnetic shock wave theory of super- 
nova origin is a definite prediction regarding an energetic electromag- 
netic pulse being emitted from the surface layer during this layer's 
ejection and subsequent adiabatic expansion. The detection of such a 
pulse and identification of it with an optically observed supernova is 
an important test of this theory of the origin of cosmic rays. 

An alternative— and less well defined— picture of cosmic ray ac- 
celeration in a supernova is that the cosmic rays are accelerated by 
the moving magnetic fields in the general turbulence following the 
explosion. General arguments can be made to explain how the energy 
spectrum might be obtained. There remains the problem of explaining 
how the highest energies could be obtained in this case, and there are 
also other detailed aspects, such as the electron acceleration and de- 
gree of energy loss, which have not been satisfactorily resolved. No 
intense, high energy gamma ray pulse similar to that of the shock 
wave acceleration theory is predicted in this case. 


PROPAGATION 

The cosmic ray energy spectra which are observed at the earth 
represent the source spectra after they have passed through inter- 
stellar matter and have been modulated within the solar system. 
Whereas presumably the solar system modulation is primarily the 
result of electromagnetic fields, interstellar space is believed to 
contain enough material along the path of the particle to change ap- 
preciably the particle energy as well as the intensity of the radiation. 

In the interstellar medium, it is normally assumed that the intensity 
in an energy interval is changed significantly only by fragmentation 
in interactions and by ionization energy loss, and not by acceleration 
nor by the complicated time-dependent magnetic effects which probably 
cause the intensity variation in the solar system. 
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If the acceleration in interstellar space is negligible, it is pos- 
sible to calculate the energy dependence of the relative abundances 
of various particle groups outside the solar system, assuming various 
source spectra, provided that the collision cross sections in inter- 
stellar space and the amount of material traversed are sufficiently 
well known. The exact nature of the solar modulation is not yet known, 
but the general belief is that it probably depends only on the velocity 
and charge-to-mass ratio of the particle. Therefore, although nuclei 
with the same charge-to-mass ratio but different charges will lose 
energy at different rates in interstellar space, the fluxes of these 
particles will be modulated in the same way, thereby permitting the 
separation of modulation effects from interstellar energy loss and 
fragmentation effects. 

In making a calculation, a particular set of simplifying assumptions 
is often made. These are: 

1. The source energy per nucleon spectra of all multiply charged 
nuclei have the same shape, at least above 100 MeV/nucleon. 
(Note that this effectively permits the spectra to be both ve- 
locity and rigidity dependent, since all of the multiply charged 
nuclei of interest at the source have nearly the same charge- 
to-mass ratio.) 

2. The abundances of both He 3 and light nuclei (3 < Z < 5) at the 
source are negligible compared to He 4 and medium nuclei 
respectively. 

3. The average interstellar mean free path is independent of the 
energy per nucleon of the particles. 

Before outlining the general procedure, these assumptions will be 
discussed. 

The similarity in source spectral shape, at least for particles of 
the same charge-to-mass ratio, is suggested by the predictions of the 
basic acceleration mechanisms, such as the simple Fermi theory, 
the principle of equipartition of energy, and the shock-wave theory 
of acceleration by supernovae. Some of the more complicated varia- 
tions of these theories can lead to predictions of differences in the en- 
ergy spectra of different charge-to-mass ratios, but this point will not 
be pursued here. Further support to the choice of spectra of the same 
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shape for particles of the same charge-to-mass ratio is obtained from 
the study of solar particles. 

With respect to the composition, the only important assumption 
is that light nuclei and He 3 are essentially absent in the source. For 
the light nuclei, this assumption is based on the fact that light nuclei 
are very rare in the universe (about 10 3 of the abundance of medium 
nuclei) because they are unstable at the high temperatures of most 
stars. The justification for assuming there is virtually no He 3 at the 
source is similar. 

Finally, the assumption that the average mean free path is in- 
dependent of energy is based on simplicity and the fact that in a simple 
steady state situation, the path length does not vary with velocity. A 
rigidity dependent average path length is possible if, for example, 
high rigidity particles escape more easily; and this alternative pos- 
sibility has been considered in the literature. It is sometimes as- 
sumed that the cosmic radiation may pass through part of the ma- 
terial in the source region itself. However, only relatively minor 
differences result from the different natures of the material in the 
source and in interstellar space. It is expected that there will be a 
distribution in path lengths about an average path length. Again, for 
most path length distributions, the results to be considered here do 
not differ much from the results of assuming that all cosmic rays 
pass through the same amount of material. 

As has already been mentioned, the material in interstellar space 
affects the cosmic ray particles by causing energy loss and fragmenta- 
tion. In order to take both of these phenomena into account, it is 
necessary to begin with the appropriate transport equation. The 
fundamental equation can be written in the form 

^ [w. (E) j. (E, x)] - w. (E)H. (E, x) , (1) 


where xis the position along the particle path; w. = (dE/dx) i ; i re- 
fers to the particular nuclear species; j i is the differential directional 
intensity per unit energy per nucleon E, for particles of type i ; and 
H. is the number of particles added or subtracted per unit volume, 
time, solid angle, and energy per nucleon. 
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In the case under consideration H. consists of two parts, the 
source term S. and the loss term due to collisions. Hence, 

H i (E, x) = S. (E, x) - j. (E, x)/A. (E) . (2) 

Here, A. is the loss mean free path, which is related to the interaction 
mean free path A . by the equation 

!/Ai = (l-Pii)^i. (3) 

where P u is the average number of particles of type i formed in the 
interaction of a type i nucleus. In turn, is given by 


S i (E, x) = 2_] i k ( E > X ) /A ki ( £ ) - (4) 

k>i 

where A k . is the mean free path for production of i-type particles 
from k-type particles. The sum in (4) can be restricted to particles 
heavier than i, because the lighter ones will not contribute to i nuclei 
in an interaction, and the i-to-i type interaction is already included 
in the last term of (2). 

Substituting (2) into (1) and rearranging terms yields 

-£ [w L j J + w. j. (1/A.) = w. S. . (5) 


Multiplying (5) by the integrating factor exp (x/A) and rewriting the 
left-hand side yields 

cE [exp (x/A.) j. wj - exp(x/A i ) w. S i . (6) 

Then (6) may be used to propagate the particle energy spectra at the 
source, j . (E, 0), through interstellar matter in small steps. The 
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step size must be sufficiently small that the variations in terms such 
as w. (E) and A. (E) introduce only a negligible error into the calculation. 

The results of such a calculation show that it is not possible to 
obtain full agreement between the calculated abundance ratios of vari- 
ous nuclear species as a function of energy/nucleon and the observed 
values. Possible explanations which have been suggested recently 
include a very strong rigidity effect in the solar modulation; a rigidity 
dependent path length; fragmentation at the source before most of the 
acceleration; and two-source models. Only the last seems not to have 
strong objections. One suggested two-source model includes a con- 
tribution at low energies from celestial objects less active than 
supernovae, and another suggests a recent close supernova burst to 
be added to the general cosmic radiation. At present the problem is 
not finally resolved, but it still appears that the interstellar medium 
probably plays the general role described above and also is responsible 
for the light nuclei and other presumed secondary products observed 
in the cosmic rays. 
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III. COSMIC-RAY ELECTRONS AND RELATED 
ASTROPHYSICAL PROBLEMS* 

R. Ramaty 
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Greenbelt, Maryland 


INTRODUCTION 

Various quantities which are relevant to our under standing of 
cosmic-ray electrons can be measured at the earth. Among these 
are 

a. The energy spectrum and intensity of the electron flux (the 
latter quantity is measured in such units as electrons/m 2 -sec- 
sr-BeV); 

b. The charge composition, i.e., the relative numbers of elec- 
trons and positrons; and 

c. Time variations. 

The effects of solar modulation, so evident in the proton spectrum, 
have not yet been detected in the electron component. However, the 
upper limits which have been established are consistent with some of 
the theoretical predictions. Other measurements which are indirectly 
useful in understanding the electron component include the nonthermal 
radio background and the spectrum of high-energy photons (X and 
gamma radiation). It is hoped that such measurements will further 
our understanding of some of the outstanding problems in high-energy 
astrophysics. Some of these problems are discussed below and in- 
clude the origin of cosmic rays, solar modulation, galactic propaga- 
tion, the connection with electromagnetic radiation, and the problem 
of energy loss, which is discussed first. 

Since the charge-to-mass ratio of electrons is so much greater 
than that of protons and other heavy nuclei, the electrons lose energy 
at a proportionately higher rate. This energy loss results in X, 
gamma, and radio radiation. The primary interactions which give 


‘Notes taken by M. L. Goldstein. 
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rise to these photons are the synchrotron, bremsstrahlung, and in- 
verse Compton processes. 


SYNCHROTRON RADIATION 


Consider an electron moving in a static , uniform magnetic field B. 
The pitch angle of the electron is 4>, and the angle made by the 
field direction and the observation direction is 8 (Figure 1). 


The power radiated by such an electron 
into solid angle dfl per unit frequency interval 
dans 



Figure 1 -Relationship of 
electron pitch angle and 
angle made by field di- 
rection and observer 
direction. 
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This expression is independent of azimuth. In Equation 1, 


SO) 

P = -C— +/3 CUCOS0 , (2) 

S / ^ 

y = E/mc 2 , E being the electron energy , (3) 


/3 is the magnitude of the component of electron velocity (* c) which 
is 1 orthogonal to B, is the component of electron velocity (* c) 
which is parallel to B, 


CO 


B 


eB/mc , 


(4) 
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and 


s/3 sin 9 sin </> 

1-/3 cos 9 cos 4> ' (5) 


When the electron becomes relativistic (y » 1), one can show that no 
appreciable power is emitted by the electron unless (&-<!>) - 1/y. 

If one were to integrate 


d 2 P 

dfi d co 


over all solid angles, the result would be equivalent to the emission 
from a region of tangled magnetic field. This integral can be done 
for electrons moving in circular orbits to yield 


-TL?JL. , n _ dP 

dfi dv dn - dv ’ (6) 


where 


dP _ e 4 B 2 J. V-i (v\ ( Sv B \ 

dv ~ m 2 c 3 y Lu 7 \ v e) T ~~y~J ’ (7) 

s 7 


Jv\ 4rr 1 

f 1 

= ilyfi 

S/3 2 J 2 ' s (2s/3)- s 2 (1-/3 2 ) J 2s (2sx) dx 

•'O 
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and 


V 


c 



(9) 


F(v/v c ) is plotted as a function of v/v c in Figure 2 for electrons of 
various values of 7. 


As y increases, the harmonics lie progressively closer together. 
The curve 7 = ro is commonly used in synchrotron calculations. It 
can be shown that as 7 “ “ , 





K s / 3 (y)dy 



Figure 2-F(v/v e ) versus v/v c for various values of J. 
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The peak in F(v/v c ) is at v/v Q % 0.3. If there is a plasma present 
with dielectric constant n 2 3 (l - w^/w 2 ) 1 /2 , where is the plasma 
frequency (4-n-Ne 2 /m ) 1/2 , and N is the electron number density, the 
curves of ¥(v/v c ) bend over at low frequencies. The effect of the 
ambient plasma, therefore, leads to a low-frequency cutoff of syn- 
chrotron radiation. The total power radiated by an electron in a 
vacuum is obtained by integrating dP/dv over all frequencies and is 
given by 


J^dz, = E = |cr 0 2 B 2 ^ 2 y 2 , (10) 

where r Q is the classical radius of the electron, e 2 /mc 2 . 

INVERSE COMPTON EFFECT 

e + y - e + y 

high-energy low-energy scattered high-energy 

electron photon electron photon 

If we let e be the energy of the incident low-energy photon and e ' 
be that of the scattered photon, then neglecting angular terms, one 
can easily show from the kinematics of the interaction that 


The rate of energy loss is just the energy lost per collision multiplied 
by the number of collisions per second. Thus, the total radiated power 
is given by 



E 


for 


7 £ 5 x 10 4 ,(11) 
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where the restriction in y values requires that classical Thomson 
scattering occurs so that ye/ me 2 « 1, and 


cr 


c 


8 

3 77 T 


2 

0 ' 


Thus synchrotron and Compton energy losses are both proportional 
to the energy density of the ambient fields. This result can be de- 
rived more directly and exactly by starting with the formula for in- 
stantaneous power radiated by an electron in an arbitrary electro- 
magnetic field (References 1 and 2), which is 


P 


2e 4 (E + Bx H) 2 - (B - E) 2 
3m 2 c 3 1 — /3 2 


( 12 ) 


Using the fact that for any energy flow which is isotropic, 
\fi ~ (Ex H)) = 0, and that for unpolarized fields 

<(/3xE) 2 > av = f/3 2 , <(/?x H) 2 ) av = |/3 2 H 2 , 


we have 



In the situation in which the only E fields present are radiation fields, 
the energy flux incident on the electron is (8/3) -n r Q 2 (cE 2 /47t). Thus 
the average rate of change of energy is 


3. 555 • • • v r Q 2 c 


(E 2 + H 2 ) 

877 - 



(de/dt> AV 


(14) 
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where E 2 /87 t is just " ph in the Compton process, and (P/mc 2 ) 2 = y 2 ,< 5 2 _ 
However, in this derivation, we have included static magnetic fields as 
well and therefore have included synchrotron losses. 


PRODUCTION OF POSITRONS AND NEGATONS IN COSMIC-RAY 
INTERACTIONS IN INTERSTELLAR SPACE 

Cosmic-ray protons undergo various interactions with inter- 
stellar hydrogen. Among these interactions are 


p + p -* p + p + a[rr + + v ) + bi7° 

p + n + rr + + a(rr* + 77 “) + b' 


2n + 2v + + a (rr + + 77 j + b77° 


The threshold for pion production in the lab system is 287 MeV. As 
the incident energy is increased above threshold, the first reaction to 
occur is pp -* pn 7 T + . At low energies, therefore, the ratio of 77 + to 77“ 
is greater than unity. At high energies, this ratio tends toward unity. 
If we call q w (T w ) the rate of production of pions per unit volume per 
unit energy of the pion (T w ), then 


%{\) = 477 JdTp j ( T p ) a-(Tp) Nf ( T p , T^) . (15) 


The units of the right-hand side are, reading from left to right, (4'") 

= steradians, [j (T p ) ] = protons/cm 2 -sec-sr-energy interval, ( T p )] 

= cm 2 , [N] = number of hydrogen atoms per cm 3 , [f (T p , T w )] 

= proton energy/pion energy. Cosmic-ray experiments directly 
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measure j (T ) , while N is derived from 21-cm radio data; f and 
cr (T ) are measured in the laboratory. 

We are interested in pion production because charged pions de- 
cay first into muons and finally into electrons via 


77 1 - /X* + V 

jjr — e ± + v + v . 


The decay energies of the two-body decay are uniquely determined 
from energy-momentum conservation although this is not true of the 
three-body decay, and thus the electrons and neutrinos are produced 
with a continuum of energies. If we let P {y * > cos &*) dy* d cos 9* 
be the probability that an electron is produced with energy y e in the 
direction 9* (where * indicates the muon's rest system), then 


q (y^) dy^ P(y e * - cos e *) dy e * d cos e * - dy e dy e * dy^q (y e , y* , y^ , (16) 

where q ly ^ ) is the number of muons of energy y p _ produced per cm 3 
per second per unit y , and 


y e - y M v* + /y A 2 - i /y e * 2 - i cos e* 


(17) 


Since we are not interested in y* and y , we define 


Q(y e ) 


jjq(y e ; y e *. yj dy e * d 7fl 

J dy e * dy^ q M (yj P (y e * , cos 9*) 


d (cos 9*) 
dy e 


(18) 
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where 


d (cos 0*) 

dy 

' e 


1 

V V - 1 +y* 2 - 1 


(19) 


One assumes that pion decay is not very sensitive to angles. Be- 
cause the 77, // mass difference is small, one can show that q., (y ;j 1 
% fin)- Analogous simplifications are difficult to make for muon 
decay because the u , e mass difference is large (--200), and the decay 
is sensitive to the angle of emission of the electrons, positrons, and 
neutrinos. The equilibrium flux of electrons in interstellar space can 
be obtained by solving a diffusion-energy loss equation. This is dis- 
cussed later. 

The measurements of the spectrum and charge composition of 
cosmic-ray electrons at the earth are shown in Figure 3. The curves 
of e s and e s + correspond to the calculated fluxes in interstellar space 
of total secondary electrons and positrons respectively. The curves 
e s + (Mod !) and e s + (Mod II) represent the modulated positron fluxes 
at the earth, using either one of the mbdulating functions given in the 
figure. As can be seen, the existing positron data are consistent with 
both models, but measurements of positrons below 100 MeV could 
easily determine the amount of electron modulation at these energies. 

The modulation of electrons and their extrapolated interstellar 
spectrum were discussed by Ramaty and Lingenfelter (Reference 3). 


PROPAGATION OF ELECTRONS IN INTERSTELLAR SPACE 

Let Q(E, r, t) be the number of electrons produced at rwith 
energy E per unit volume per second. While losing energy, particles 
diffuse away from this source region. The equation which describes 
this is 


5u 

It 


- V 


(DVu) +-^Y (Eu) - 


Q(E, r) , 


( 20 ) 
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ELECTRON ENERGY (BeV) 


Figure 3-Earth measurements of the spectrum and charg- 
composition of cosmic-ray electrons. 
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where u is the electron density, Dthe diffusion coefficient, and E(E) 
is the rate of energy loss caused by ionization, bremsstrahlung, 
synchrotron, and inverse Compton losses. Furthermore, we restrict 
ourselves to the example of an infinite volume. The solution is 


U (E, r) 



- | r-r ' | 2 /4 Dt 

e 

(4ttDt) 3/2 


Q(E', r') 


( 21 ) 



( 22 ) 


We assume 


Q(E, r) - q(E)p(r) , 


(23) 


so that 


u(E, r) 


JdE' 


q(E') {(£, r) 


(24) 


where 



- | r- r ' | 2 Adt 


3/2 


r) 


(4v Dr) 


P(r') , 


(25) 
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and 


a 2 /2D 


where a is the radius of the diffusing region (e.g., the galaxy). 

A simplified result which at least qualitatively illustrates some 
features of the more exact solutions may be obtained by setting 


V • (DVu) = u/T . 


In this case, f (£) = e f , and Tis the time required for a cosmic ray 
to diffuse out of the galaxy. If p is the density of matter in interstellar 
space measured in gm/cm 2 , we can define a path length x - pcT. 
Synchrotron and Compton energy losses can be combined by setting 


dy 

dt 


KJ 


2 


The characteristic time for energy loss is then 



The solution of the diffusion equation now tells us that after a 
time T all electrons having energies greater than y 0 , where *7 0 - 1/T, 
will have lost a significant fraction of their energy. For the special 
case in which y" n , the solution is illustrated in Figure 4. 

The fact that no break in the measured electron spectrum is seen 
up to about 300 BeV indicates that T < 10 6 years. Combining this in- 
formation with estimates of the path length X determined from studies 
of the isotopes H 2 and He 3 as well as the nuclei Li, Be, and B, one 
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LOG 7 


Figure 4 — Solution to special case Q a r n . 


finds that p % lg/cm 2 . This indicates that the high-energy electrons 
must be confined in the galactic disk. Parker (Reference 4) arrived 
at a similar conclusion for the majority of cosmic rays from dynami- 
cal considerations of the cosmic rays and the interstellar gas and 
magnetic fields. 

Other questions of astrophysical interest related to cosmic- ray 
electrons are the origin of the nonthermal radio background and the 
diffuse components of X- and gamma rays. These questions were not 
discussed in this lecture and the reader is referred to the references 
cited. 
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IV. COSMIC-RAY PHOTONS 


E. Boldt 
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THERMAL SOURCES 

Unlike laboratory situations in physics observations, photon astro- 
physics is limited to remote observations without the experimental 
capability of controlling and perturbing the physical conditions being in- 
vestigated. Hence, we shall attempt to interpret the production and 
propagation of the photons of astronomy in terms of processes already 
examined in the laboratory (e.g., synchrontron, bremsstrahlung, and 
black-body radiation; Compton scattering, photoelectric absorption, 
etc.). 


The source of the most intense flux of cosmic photons observed at 
the earth is the sun. For wavelengths from~10 3 A to 1 cm, this source 
appears as a black body of T a* 6 x 10 3 °K. The spectrum peaks in the 
visible region, where the earth's atmosphere is relatively transparent. 
For radio waves exceeding ~1 cm wavelength, the intensity exceeds that 
of a black body at 6 x 10 3 °K and approaches that of a black body at 10 6o K 
for very long wavelengths (greater than^lm). The emission in the ul- 
traviolet ($10 3 A) and X-ray regions (<10 2 A) also exceeds that of a 
black body at 6 x 10 3o K but does not, in fact, follow a black -body spec- 
trum at an elevated temperature. The X-ray emission from the quiet 
sun indicates a hot (>10 6 °K), transparent thermal electron plasma; 
solar flares emit hard X-rays ($10 keV). In no instances have solar 
gamma rays ($MeV photons) yet been observed. 

Since thermal radiation appears to be an important astrophysical 
phenomenon, it is useful to relate observations to the black-body model. 
The number of photons N per degree of freedom within a volume F of 
phase space is given by 


N 


f f ( p , x) d 3 p d 3 x 

J h3 ’ 


(1) 
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where p and x are the photon momentum and position vectors, respec- 
tively, and h is Planck’s constant. The function f ( p , x) is the distri- 
bution function for the photons of a disordered (e.g., thermal) radiation 
field, obeying Bose statistics. There are two degrees of freedom, one 
for each of the two spin states available to a photon. If we consider a 
small region T in the vicinity of (p, x) and measure N in several tem- 
porally separate, independent observations, we obtain 

J ) 

— - 1 + f(p> x) . (2) 

From Equation 2 we note that thermal photon fluctuations are some- 
what similar to what would be expected from the classical statistics 
associated with particle counting but are in fact always greater by an 
amount given by f . This dual behavior for the statistics of photons, 
whereby they behave somewhat as particles on one hand, but yet ex- 
hibit the "bunching" characteristic of wave interference on the other 
hand, was first noted by Einstein for the black-body distribution func- 
tion, viz: 



^ b £exp (hv/kT) - lj 


( 3 ) 


For this black-body distribution function, the photon bunching (Equa- 
tion 2). exhibits two interesting limits as follows: 


((n-< n >) 2 ) 

(for | 

[a 

* ’) ■ 

(4a) 

((»-<">)’) kT 

(for 

r hv i 

LkT. 

I* 1 )- 

(4b) 


Table 1 summarizes the amount of photon bunching at various 
characteristic wavelengths for several important astronomical sources 
ranging from the 3°K universal microwave background radiation to the 
10 6 °K radio emission of the sun. For this range, there is appreciable 
bunching for A. > 1 cm and almost none at shorter wavelengths (e.g., 



Table 1. 


Photon Bunching, 


((n- <N>) 2 y<N> 


Radiation field 

kT 

Wavelength 

(eV) 

5 x 10 3 A 

10 2 cm 

1 cm 

21 cm 

1 m 

Universal microwave 
(3°K) 

3 x 10 -4 


(1) 

(3) 

(40) 


Interstellar hydrogen 
(100°K) (21 cm) 

9 x 10 - 3 




(10 3 ) 


Optical sun 
(6 x 10 3 o K) 

5 x 10 - 1 

(1.01) 


(4 x 103) 



Optical B-star 
(14 x 10 3 o K) 

1 

(1.2) 


(104) 



Radio sun 
(10 6 °K) 

90 





(9 x 107) 

Photon energy 

3 eV 

10~ 2 eV 

10~ 4 eV 

6 x i(T 6 eV 

10' 6 eV 


Ol 
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optical). Thus, there is a rather clear separation here between (1) 
radio astronomy as a carrier of wave information about the radiation 
field observed, and (2) optical and X-ray astronomy as a carrier of 
particle-type information, in the sense that each observed quantum is 
statistically independent of all others. For a nonequilibrium radiation 
field, the photon bunching may be very different from that for a thermal 
field. As an example of a nonequilibrium field, consider an ideal 
laser and note that in this extreme case we recover classical particle 
behavior in the sense that 


(("-<">)’) . , , 
<»> 


( 5 ) 


even though the number of photons within a single cell (h 3 ) of phase 
space could be >10 10 . 


From Equation 1, we may construct the photon intensity of a ra- 
diation field as 


dj _ 2cp 2 f ' , . 

— : — | - r — photons/(area-time-solid angle-unit p ) . 

d|p| x ' 


The spectral intensity is then obtained as 


- /u n d J 

- (hv) ^ 


2hv 3 { 


( 7 ) 


For a black body (Equation 3), the spectral intensity (Equation 7) 
becomes 


I 


V 



exp (+hv/kT) 


- 1 ] 


(8a) 



(as , 


I 


(8b) 
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( hv 

V as t?r - 0 


(8c) 


Equation 8c is the Rayleigh-Jeans limit. The "brightness" tem- 
perature (T b ) to be associated with an observed spectral intensity I v 
is given via Equation 8 as 


kT B 


hv 


- tn 


2 hv 3 



(all h^) , 


(9a) 


I 

kT B % -i \ 2 



hv 3 



(9b) 


If in fact the observed source is thermal and the observations are 
in the regime (kT/hv) » 1 , then Equation 9b is an appropriate approxi- 
mation, the one usually used in radio astronomy, and there should be 
an appreciable bunching of photons given via Equation 4b as 

(( N -< N >) 2 ) kT c 2 I v 

(N) hv 2hv 3 " ( 10 > 

In the other extreme, if one inadvertently attributes a brightness tem- 
perature to a laser source, then the lack of photon bunching would in- 
dicate a nonequilibrium aspect of the source. Thus, for the radio 
astronomy regime ( c 2 l v » hv 3 ) , we can identify a thermal source by 
observing over a narrow band at v and comparing the measures of 
photon bunching and I v with the relation fixed by Equation 10. For the 
optical and X-ray regime, this thermal identification at a spectral 
point is unavailable in most instances. 

For a black-body source, the measured brightness temperature is 
of course independent of the observed wavelength. For the observed 
extended sources of radio-emission (e.g., Galaxy, M31, Cas A, Cent A), 
the brightness temperature generally decreases as the observed wave- 
length decreases. Two notable exceptions to this are (1) the universal 
microwave background, which exhibits an apparently constant brightness 
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temperature of ~3°K over almost two decades of wavelength, and (2) 
the sun, which gives a rather constant 10 6o K for wavelengths longer 
than 1 meter. 

What is the brightness temperature of the Crab Nebula in the 
X-ray region? To evaluate this, we use the recently obtained X-ray 
information that 


10 11 eV 

v ( Hz - cm 2 -sec-sr) 


( 11 ) 


over a solid angle ^10 7 sr. 

Using Equation 11 in Equation 9a, wte obtain for the Crab Nebula 


hv o 

= % 60 (aU% 1A) . (12) 


If the Crab is a thermal source of X-rays then, from Equation 12, 
we note that it is definitely in the particle (i.e., nonbunching photons) 
regime at a high brightness temperature (T £10 6 °K). As we shall dis- 
cuss later, the X-ray spectrum of the Crab (Equation 11) is not of the 
form expected for an isothermal source, be it a black body or a trans- 
parent, hot plasma radiating via the bremsstrahlung mechanism. 

The sun and Sco X-l, the brightest X-ray star yet observed, both 
exhibit X-ray spectra consistent with that expected from hot thermal 
plasmas at £l0 6o K and^l0 7o K, respectively. The photon distribution 
function f inferred from Equation 7 is definitely much less than unity 
for the solar X-rays, indicating that even for a clearly thermal source 
the X-rays are again in the particle regime; the corresponding analysis 
for Sco X-l still lacks the information required on the lower limit to 
the angular size. 

For a black body, the isotropic intensity I of energy is 



where a = Stefan-Boltzmann constant. 
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Table 2 shows the total output for several objects of interest. For 
black-body radiation, Equation 13 was used to evaluate ( 2nr ) 2 I for the 
total optical output of the sun, with r as the solar radius, and for the 
total microwave power exchanged by the galaxy over the spherical halo 
(r ~10 23 cm). Although at only 3°K, the galactic microwave power 
equals the optical power output for ~10 n suns at a black-body tem- 
perature of 6 x lo 3 °K. The sun and Sco X-l are probably both thermal 
X-ray sources, but the output of Sco X-l equals that of ~10 13 suns as 
far as X-rays are concerned. These comparisons illustrate the wide 
range of thermal situations that obtain in astrophysics. 


Table 2. Source Strengths 


Source 

Output 
(ergs/ sec) 

X-ray sun 

10 23 

(Thermal plasma £10 6 o K) 


Optical sun 

4 x lo 33 

(Black body, 6 x 10 3 °K) 


Sco X-l 

^10 33 

(Optical, blue star) 


Sco X-l 

~10 36 

(X-rays, thermal plasma, £l0 7o K) 


Crab Nebula 

-^20 36 

(Radio and optical, synchrotron radiation) 


Crab Nebula 

~10 37 

(X-rays) 


Galactic microwave 

6 x 10 44 

(3°K black body) 



The observed apparent brightness temperature of the Crab 
(T b £10 6 °K for X-rays) may be used to calculate I, via Equation 13. 
Comparing the integrated X-ray flux from the Crab (see Table 2) with 
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( 2ttt) 2 I fixes the radius r to be associated with a black-body source 
as r ^10 8 cm (i.e., smaller than the earth) whereas X-ray observa- 
tions clearly indicate a size comparable to a light year. This rules 
out black-body emission for the bulk of the X-rays from the Crab 
Nebula. 

The black-body photons of the microwave background and starlight 
constitute scatterers for high-energy charged particles, as we shall 
describe later. The scattered photons may then become energetic 
enough to be called X-rays. We shall find it useful to relate the en- 
ergy of the photons of the black-body target with the temperature by 
noting that 


(hv) = 2.7 kT (14) 

and, from Wien's Displacement Law, that the maximum T occurs at 


M. = 2.8 kT . 


(15) 


It is also important, in considering the black body as a target, to 
specify the spatial number density n of photons, given by 


n 


477 I 


20 T 3 , 


(16) 


where n is in cm -3 , and T is in °K. Values of photon density for the 
microwave background, visible light, and X-rays in various astronom- 
ical regions of interest are listed in Table 3. 


Table 3. Photon Density (cm -2 ) 


Radiation field 

Galaxy 

Crab Nebula 

Sun 

Microwave background 
(3°K black body) 

5 x 10 2 

5 x io 2 

5 x io 2 

Optical 

1 

10 2 

4 x loi 2 

X-rays 

10' 8 

10" 1 

1 
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NONEQUILIBRIUM RADIATION PROCESSES 

The photon source function q at hv for radiating charged particles 
in space is 


q(photons/cm 3 -sec-erg ) 



where W is the kinetic energy of the radiating particle, and N is the 
number of radiating particles. 

The observed photon flux is 


j( phot on s/cm 2 - sec- sr-erg ) - 4 ^ q dr , (18) 

Jo 

where r is the radial distance from the observer in the direction of the 
observed flux. 

The central problem then is to define the charged-particle popula- 
tion and the radiation losses. The interaction of an electron with its 
electromagnetic environment causes radiation loss as follows: 


” dW - /u \ 

~St - va o n v\>W - 


(19) 


where ( hv ) is the mean energy of radiated photons, v is the electron ve 
locity, a Q is the Thompson cross-section, and is the number density 
of effective target quanta in the ambient electromagnetic environment. 

The starting point of this analysis is the relativistic classical ra- 
diation law (e.g., Reference 1, Equation 9-72), viz: 


-dW 

dt 


— r 2 
3 r e 


(e + ^ xH ) 2 - ~ (E' v) 2 ] 


y 2 c 


(20) 


where r e is the classical electron radius, and y~ 2 - 1 - (v/c) 2 . 
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An environment of plane electromagnetic waves corresponds to 


H - n x E , 

(21a) 

n • E = 0 , 

(21b) 

E 2 = H 2 , 

(21c) 


where n is a unit vector along E x H . 

We consider the case where the vectors (n) are randomly oriented, 
for which we have 


( n • v ) - 0 . (22) 

This prescription (Equations 21 and 22) allows an evaluation of 
Equation 20 as 

/-dW\ _ / fi 2 \ „ 

[ dt j EH - CT 0 c O>EH Y + Tjy - (23) 

where 

H 2 + E 2 E^ _ H 2 

877 477 ~ 477 

In the relativistic limit, we have 


/-dW\ 4 

l^rj EH - 3 °- 0 c (^)eh> • (24) 

For a radiation field environment consisting of quanta of average 
energy e = (p/ n v ) , a comparison of Equation 24 with Equation 19 shows 
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that, in the relativistic limit, 


(hu ) - 3 r 2 e . (25) 

This relation, between the mean energy of the radiated photon and 
that of the quanta of the medium, describes what is known as inverse 
Compton scattering. The soft quanta of the environment are typically 
optical starlight and radio waves; these serve as targets for high- 
energy electrons. 

Another case of great astrophysical importance corresponds to a 
finite H field and no E field (i.e., E = 0); this is the condition for mag- 
netic bremsstrahlung ("synchrontron radiation"). We obtain, from 
Equation 20, 


(~ar) H - i r e 2 c/3 2 y 2 h 2 

where 


(-,)■ ■ ^ - 


For the situation of random field orientations, we use 


We note that 



(27) 


P H ~ H 2 / $7T . 


(28) 
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Then, from Equations 26, 27, and 28, we obtain 

(“ar) H = t°-() c <»h y 2 p 2 - (29a) 

[ -dW \ 4 

ldrJ H - 3 CT o c (^ H r 2 as ,5 - 1 . (29b) 


It is important to recognize that this radiation loss formula (Equa- 
tion 29b) for synchrotron emission is identical to that for inverse Comp- 
ton losses (Equation 24) in the ultrarelativistic regime, when we ex- 
change (p) H and (p)^. 

From Equations 19 and 29a, we identify 

(hv) = ^y 2 /3e , (30) 


where e = (p) H )n y . 

In the nonrelativistic limit, for cyclotron radiation, the radiated 
photon energy is 


(hv ) - hv - b ^ H x . 


(31) 


Using Equations 31 and 30 for y -1 prescribes e, as 


3 \ h me h jJ (32} 

€ - Tf £ •• 

Inserting Equation 32 into Equation 30 yields a general expression 
for the mean radiated photon energy as 


(hv) = y 2 h Jh, . 


(33) 
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For the ultrarelativistic situation, the photon energy at the peak 
of intensity is~0.3 (hi, ) (see Reference 2). 

For an electric field alone (H = 0), we have the condition for 
bremsstrahlung radiation, and Equation 20 gives 


2cr o c( »i 


(>-*)■ 


where we have again considered a randomly oriented field for which 


(cos 2 #) = 3 


where 9 is the angle defined by E and v. 

We consider that the electric field arises from a number density 
n o °f discrete charges Ze that may be localized to within an impact 
parameter b. Therefore, we construct the electrostatic energy density 
P E required for Equation 34 as follows: 


•••(c 


( Ze/r 2 ) : 


47tt 2 dr 


The smallest impact parameter b, for localization of the charge, 
is given by wave mechanics according to the prescription 


b = \/4 , 


where X - h/(mv). Insertion of Equation 36 into Equation 35 yields 


E 77 


mc 2 /3n , 
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where a ® e 2 /fTc = 1/137. Using p E in the nonrelativistic limit (0-> 0, 
y - 1) of Equation 34 yields 

(if)_ = T aZ2 r e 2 ( m 0 c2 )^ cn 0 • (38) 

E NR 

This is the same as the nonrelativistic (NR) expression given by Heitler 
(Reference 3) ; the extreme-relativistic (ER) expression is given by 
Heitler as 

[~w) = 4aZ 2 r/ (m 0 c 2 )ycn 0 [^n (l83Z‘ 1/3 ) + yg] • (39) 

E ER 

The conventional approximation to (hy) is 

(hv) = W . (40) 

Comparing Equation 38 with Equation 19, using the identification 
(Equation 40) for (hy ) , gives 


8a n o 


(41) 


Table 4 summarizes the relevant parameters for electron radiation 
losses by synchrotron emission, inverse Compton scattering, and 
bremsstrahlung radiation, as described earlier. The numerical values 
for these parameters for two important extended regions, the Crab 
Nebula and the galaxy, are listed in Table 5. The equivalent kinetic 
temperature T k and black-body temperature T fi are defined also in 
Table 5 for each of the electromagnetic states described for the en- 
vironment. For all the astronomical situations enumerated in Table 5, 
the equivalent black-body temperature is always the order of 10°K, 
whereas the equivalent kinetic temperature varies by more than 17 
orders of magnitude. The electromagnetic state listed as "Coulomb" 
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Table 4. Review of Some Radiation Mechanisms 
^Radiation Loss Formulation: - (dW/dt) = va Q (p/e) (by) .J 


Parameter 

Synchrotron 

(ultrarelativistic) 

Inverse Compton 
(ultrarelativistic ; 
for black body) 

Bremsstrahlung 

(nonrelativistic) 

(hv) 

4 2 

3 y e 

4 2 

3 y 2 e 

1 

2 W 

e 

3 e 

4 ” me 

2.7 kT 


p 

H 2 /8tt 

4ctT 4 / c 

K)^o , 

77 K c ) 

-( dW/dt) 

A - W 

4 ( P H \ 

4 f Peh \ 

3 "» cr U^j 

4 (aZ 2 ) 

TT a 0 C ~T~ n 0 


refers to the electrostatic fields generated by the indicated astronom- 
ical proton gas, which constitutes the bremsstrahlung target. 

The decay constants A for the various radiation-loss mechanisms 
that occur for electrons are presented in Table 4, and the numerical 
values for the corresponding lifetimes (i.e., A" 1 ) are listed for many 
situations of interest along with the value of y required in each in- 
stance to obtain radiation in the X-ray region (i.e., ~1A) . 

The favored process for radiative loss in the Crab Nebula is ob- 
tained from inspection of Table 6, which indicates that the lifetime of 
a 10 13 eV electron for synchrotron X-ray emission is only about 1 year. 
As indicated in Table 6, the lifetime is proportional to y ~ 1 . Hence, 
for a lifetime comparable to the age of the Crab Nebula (~10 3 years), 
it is necessary for an electron to have an energy -10 10 eV. Since the 
average energy of a synchrotron-radiated photon is proportional to 7 2 , 
an electron of 10 10 eV corresponds to photons in the far infrared. This 
is the region where the observed electromagnetic spectrum of the Crab 
Nebula exhibits a change in spectral index. Hence, the break in the 












Table 5. Some Astronomical Photon Gases and Fields ® 


Electromagnetic state 

Quanta 

Kinetic temp, 
(equivalent) 
T k = e/2. 7k 
(°K) 

Black body temp, 
(equivalent) 

T b = ( pc/4a) 1/4 
(°K) 

e 

(eV/photon) 

P 

(eV/cm 3 ) 

n„ = (p/e) 

(photon/cm 3 ) 

o 

X-ray photons (~1A) 






Galaxy 

10 4 

10“ s 

M 

O 

1 

vO 

4 x 10 7 

0.2 

Crab 

10 4 

10 3 

10“ 1 

4 x 10 7 

20 

Magnetic 






Galaxy (5 x 10~ 6 G) 

5 x 10~ 14 

1 

10 13 

© 

1 

o 

T— 1 

3 

Crab (1(T 3 G) 

io- 11 

4 x io 4 

10 15 

00 

l 

O 

T— 1 

40 

Optical photons (^6 x 10 3 A) 






Galaxy 

2 

1 

0.5 

10 4 

3 

Crab 

2 

10 2 

50 

10 4 

10 

Microwave background 






Universal (3°K) 

7 x 10 -4 

1 

10 3 

3 

3 

Coulomb 






Galaxy (n 0 ~ 1/cm 3 ) 

ySW/4 

10 3 /3 

4 x 10 3/W 

/3W/ (Ilk) 

~20/3 1/4 

Crab (n 0 ~ 1/cm 3 ) 

/3W/4 

10 3 /3 

4 x io 3/w 

/3W/ (Ilk) 

~2 0/3 1/4 
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Table 6. Conditions Required for X-ray (^1A) Production Mechanisms 
in Extended Sources, and Associated Lifetimes 








Inverse Compton 







oyncnx'otron 

Optical field 

Microwave field 

Diem o an cui± uug 

Extended source 

Required 

electron 

1/Ah 

Required 

electron 

1 A EH 

Required 

electron 

^Aeh 

Required 

electron 

1/ 


7 

He 

(cm -3) 

(years) 

y 

He 

(cm -3 ) 

(years) 

7 

n e 

(cm -3 ) 

(years) 

P 

He 

(cm -3 ) 

(years) 

Crab Nebula 
(~10 18 cm) 

10 7 

10" 10 

1 

10 2 

10 3 

10 7 

10 4 

10" 6 

10 7 

0.3 

10 7 

10 3 

Galaxy 
(~10 23 cm) 

10 8 

10' 21 

10 4 

10 2 

10“ 8 

10 9 

10 4 

10" 11 

10 7 

0.3 

10" 6 

10 3 


05 

07 
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spectrum appears to be consistent with the synchrotron picture, but the 
extremely short lifetime for X-ray synchrotron emission would require 
much continuous activity within the remnant, persisting ~ 10 3 years 
after the supernova. 

We consider now the spatial density n x of X-rays generated by 
electrons of density n e within a spherical region (radius R) of space, 
interacting with the environmental field quanta (density n„) of that re- 
gion. From Equations 17, 18, and 20, we obtain that at the origin 


n = (" ( hl ") J d 

•'o 

where 


( n e n .) 


Taking R ^ 10 18 cm for the Crab Nebula and R % 10 2 3 cm for the gal- 
axy, the electron concentration n e that is required for the observed dif- 
fuse X-ray density within each of these source regions may be calculated 
by using Equation 42 with the n v appropriate to the radiation interaction 
considered (see Table 5 for numerical values of n„). As shown in 
Table 5, the X-ray photon densities n x within the Crab Nebula and the 
galaxy are -'T0~ 1 /cm 3 and 10~ 9 /em 3 , respectively, these are the num- 
bers used in Equation 42 for evaluating the magnitudes of n e listed in 
Table 6. 


(hv) 


0 (\ e l v / J ^ ’ 






V \ 

v J."- n ' d V' 


By inspecting Table 6, we infer that the energy density for elec- 
trons of ''TO 13 eV in the Crab Nebula which would be responsible for 
X-ray emission by synchrotron radiation amounts to ~1 percent of the 
energy density in the magnetic field itself. The energy density for the 
electrons of ''TO 10 eV required for X-ray production by inverse Comp- 
ton scattering of the microwave background is ~10 percent of the en- 
ergy density of the magnetic field confining the electrons. However, 
such a density of ~10 10 eV electrons would synchrotron-radiate more 
photons in the radio band than observed. The electron number density 
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required for bremsstrahlung generation of the X-rays of the Crab 
Nebula is prohibitively high for such an extended source. Although it 
appears from this analysis that synchrotron radiation is likely, prob- 
ably the only unambiguous way to establish synchrotron radiation as 
the dominant X-ray-generating mechanism of the Crab Nebula is 
through a definitive measurement of the X-ray polarization. 

For the galaxy, Table 6 indicates that electrons of ~ 10 14 eV at an 
energy density of ~10“ 7 eV/cm 3 would give sufficient radiation by the 
synchrotron process; this density is higher than direct cosmic-ray 
electron measurements at lower energies would suggest. Inverse 
Compton scattering by the microwave background would require elec- 
trons of ~10 10 eV at^KT 1 eV/cm 3 — again too high. Inverse Compton 
scattering by starlight would require electrons of ~ 10 s eV at an energy 
density of ~1 eV/cm 3 — too high again. Electrons of ~ 1 0 4 eV would 
generate sufficient galactic X-rays by bremsstrahlung for an electron 
energy density of ''“10” 2 eV/cm 3 , comparable to the kinetic energy den- 
sity for the interstellar hydrogen gas. The effect of X-ray generation 
in the metagalaxy has not been considered here, but could be important. 

SPECTRA 

The spectral representations for the radiation-loss mechanisms 
of bremsstrahlung, inverse Compton scattering, and synchrotron emis- 
sion are approximated as follows: 

Bremsstrahlung: 



Inverse Compton: 



d 

dhz^ 

ldtJ EH 

© E [s(h^- |y 2 e) 

] - (44) 

Synchrotron: 




d 

d\w 

($L\ ^ 
'dt - 

(hf) h [8(h^-0.3y 2 <0 

1 • (45) 
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(For a discussion of the approximation given by Equation 45, see Ref- 
erence 2.) 

We may now construct the photon source functions for bremsstrah- 
lung, inverse Compton, and synchrotron radiation by evaluating Equa- 
tion 17 for Equations 43, 44, and 45, respectively, using the basic ra- 
diation loss formula (Equation 19) : 


where (h^) 
tion 41) ; 


where 


and 


where 



= W/2 (see Equation 40) and n v = 8a/ tt (n Q /3 2 ) (see Equa- 


( dW \ 

d 

/d 3 N\ 


IdEjj 

dW 

\dx 3 / ’ 

(47) 


7 2 ~ J (hu/e) ; 


/ dW \ d / d 3 N \ 

% - c n v VdhjjJ "3W / < 48 > 


o / 3 e 

y 2 % 3(h^/e) and e - 4; fr — H x . 
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We examine here the functional dependence of q on photon energy 
for the case of a power-law electron distribution: 



(49) 


Using Equation 49 for q as given by Equations 46, 47, and 48 gives 


q E « 

(hi/)"" , 

(50) 

^EH 

(hy) _ < n+ U/2 f 

(51) 

% “ 

(h^)-(" +1 >/2 . 

(52) 


It is interesting to note that the power law index for relativistic, 
charged cosmic rays is observed to be n % 5/2, and the diffuse cosmic 
photon spectrum exhibits a power law above 10 4 eV of index ~2. These 
observations are then spectrally consistent with the emission mecha- 
nisms of synchrotron radiation and inverse Compton scattering as in- 
dicated in Equations 51 and 52 where (n + l)/2 %2, for n Zs/2. The 
X-ray spectrum of the Crab Nebula also exhibits a power law spectrum 
of index 2. 

Another electron spectrum of astrophysical interest is that for a 
thermal gas, viz: 


d 

dW 



(W) 1/2 exp(-W/kT) , 


(53) 


where T is the kinetic temperature for the Maxwellian distribution. 
These hot electrons will radiate by bremsstrahlung collisions with 
the atoms and positive ions of the gas. If one uses the Maxwellian 
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distribution for electrons in evaluating q E , then 



where n and n„ are the number densities of electrons and ions (charge 

e U 

Ze), respectively. 

The quantity g in Equation 54 is the free-free Gaunt factor, which 
is about unity for most atoms, and which is given explicitly for hydro- 
gen by 


A 

g =■ — [exp (hv/ 2kT>] K Q (hv/2kT) , (55) 

where K 0 (x) is the modified Bessel function of second kind (Drummond, 
1961). The total power radiated per unit volume by thermal bremsstrah- 
lung is obtained by integrating q E as given by Equation 54 over all pho- 
ton energies and yields 


J n 


(hv) q dhv - 14 


10 


-24 




J-) 

10 ‘I 


1/2 


ergs 


cm sec 


(56) 


where T is in °K. 

As an example of a thermal source, consider the X-ray star Sco 
X-l, which exhibits a spectrum of the form (Equation 54) to be ex- 
pected from a transparent thermal plasma at T < 10 So K. The total out- 
put in X-rays is estimated as '-10 36 ergs/sec. If we take the volume 
as comparable to the sun (~10 33 cm 3 ), we may use Equation 56 to eval- 
uate (n e n 0 ) as ''TO 26 cm -6 . For the case n g = n Q , the number density 
of electrons and ions in this hot gas is ~10 13 cm -3 . 

Now that the thermal state of Sco X-l has been specified, we may 
test if the associated hot plasma is indeed transparent. To do this, we 
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need the absorption cross-section (a a ) , which may be obtained from 
Kirchhoff's law: 


- CCr 2 n e B v • ( 57 ) 

where B v is the number of photons per unit volume per unit photon 
energy interval for a black body at temperature T . 

From Equations 1 and 3, 


4 77 (h^) 2 

(he) 3 [exp ( hiVkT ) - l] 


( 58 ) 


We may now evaluate the absorption mean free path [a. = ( n o- a ) 1 
from Equations 54, 57, and 58 and obtain, for hi>>kT, 


X, 


4 x 10 38 
gZ 2 n e n o 



( 59 ) 


where X is in cm, hv is in eV, n is in cm 3 , and T is in °K. 

For X-ray emission f-lA) from Sco X-l, we use Equation 59 to 
determine that k £ 10 24 cm, which is larger than the galaxy. Hence, 
the plasma of Sco X-l is transparent, and our picture remains consistent. 

The dominant mechanism for X-ray absorption is the photoelectric 
effect. For light elements (Z <10), the absorption length at the energy 
of the K -edge (51 keV) is 510 18 atoms/cm 2 . For unit atom number 
density, this absorption length corresponds to 1 light year. Hence, 
photoelectric absorption by the interstellar medium causes an appre- 
ciable attenuation of stellar X-rays for energies less than 1 keV. 

COSMIC PHOTONS FROM PROTON INTERACTIONS 
X-Rays 

In ionizing the atomic hydrogen of the interstellar gas, a cosmic- 
ray proton sometimes imparts appreciable energy to the ejected 
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electron; in such a case, the electron is called a "knock-on electron." 
This electron undergoes the acceleration required for attaining its 
exit velocity during the collision and therefore radiates through the 
bremsstrahlung mechanism, which in this instance is called "inner 
bremsstrahlung" (Reference 5). For nonrelativistic collisions, the 
energy imparted to an electron by a proton of kinetic energy U is 
limited by 


W £ 4 fe) U ’ ( 60 ) 

where m Q is the electron rest mass, and M 0 is the proton rest mass. 

The source function for inner bremsstrahlung (q r ) is obtained 
from 


hv qj 



where dW/dhv gives the spectral energy radiation of the knock-on elec- 
tron as 


dW 

dim 


2 a 

3 7 T 


(*.V 


(62) 


and do-/ dW is the differential cross-section for producing a knock-on 
electron of kinetic energy W (velocity /3 e ) with a proton of kinetic en- 


477 ( r g ) 2 

(0 e ) 2 te>) Sw ' (63) 

In Equation 61 dn/dU is the spatial spectral density of cosmic-ray pro- 
tons atu, ionizing hydrogen gas of n 0 atoms per unit volume. 


ergy U (velocity 0 ) ; thus 


dcr 

dW 
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The inner integral of Equation 61 may be evaluated with Equa- 
tions 62 and 63 as 



This is similar in form to the expression for the ionization loss 
per unit path length du/dx, given by 


dU 

dx 


3 CT o( m o c2 H n 
2 n o 


2m 0 


0 2 




(65) 


where I is the ionization potential for the atom. 

Using Equation 65 in Equation 61, via Equation 64, we obtain 


(hv) q : 


( m o c2 ) 


I 


in 


/2m. c 2 0 2 


hv 


u= ( M o/ 4m o) hi/ 


2m 0 c 


dj dU 

2 p 2 \ dU dx 


dU 


, (66a) 


(h^) q t 



(66b) 


where (dj/dU) is the spectral intensity of cosmic-ray protons/(cm 2 - 
sr-sec-eV) . 


This simple form (Equation 66b) for the inner bremsstrahlung 
source function allows us to make a comparison with the ionization en- 
ergy dissipated per second per unit volume of the interstellar hydrogen 
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gas by all charged cosmic rays above energy U, defined by 


Q = 



(67) 


The Q required to keep the interstellar gas heated to the observed 
temperature of 100°K is Q ^ 5 x 10“ 14 eV/cm 3 -sec, for a hydrogen 
density of 1 atom/cm 3 . This heating is likely to arise from low-energy 
cosmic rays (U>15 MeV). (See Reference 6.) 

We may write the source function for inner bremsstrahlung in 
terms of Q as follows: 


2 a Q n m o\ 

(hl/)q i % for 1 vv u o ’ ( 68 > 


where U Q is the lowest cosmic-ray proton energy that makes an ap- 
preciable contribution toQ. Using U Q = 15 * 10 6 eV and Q = 5 x 10' 14 
eV/cm 3 -sec, we obtain from Equation 68 that inner bremsstrahlung 
generation is (h^ q t ) % 10' 22 eV/eV-cm 2 -sec, for hv <(4m 0/ /u o ) u 
(= 3 x io 4 eV). For a direction that intercepts a depth -t of inter- 
stellar gas, the corresponding inner bremsstrahlung X-ray intensity I 
would be 


I 


477 (^l) 


eV 

eV cm 2 -sec-sr 


(69) 


For i = 10 23 cm, the diameter of the galactic disk, we get I % 1 
eV/cm 2 -sec-eV. This is comparable to the observed background flux 
at hv ~ 30 keV. If charged cosmic rays of U > 15 MeV are indeed re- 
sponsible for the heating of the interstellar gas, then we should be able 
to track the interstellar gas distribution with a significant component 
of the observed diffuse X-ray flux at ~30 keV. The spectral shape of 
this X-ray emission is an integral measure (Equation 66b) of the low- 
energy, charged cosmic-ray spectrum. Because of the drastic solar 
modulation of such low-energy cosmic rays, the inner bremsstrahlung 
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X-rays generated in interstellar hydrogen might provide the only 
available means for such spectroscopy. 

Consider an interstellar charged, cosmic-ray spectrum of the 
form 



U" n . 


(70) 


Using Equations 70 and 65 in Equation 66b gives 


(Mqj ~ (hv) n . ( 71 ) 

These expressions (Equations 70 and 71) indicate the basic aspect of 
interstellar inner bremsstrahlung spectroscopy for charged cosmic 
rays at low energy. 

Gamma Rays 

At U it 1 GeV, cosmic-ray protons produce pions via nuclear inter- 
actions with the hydrogen of the interstellar medium. The neutral pion 
decays into two photons with energies in the range 


K c2 ( 1 + ^Y 2 M „c 2 /'i-zVV 72 

— {—j >hv>— -\rr^j ’ ( 72 ) 

where /3 n = velocity of the pion (mass M„) in units of c and c 2 a; 140 
MeV. The mean and most probable energy for each photon is E^/2, 
where is the total pion energy. Near pion -production threshold. 



where E p is the total cosmic-ray proton energy, and M is the proton 
rest mass. 
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For an ultrarelativistic case, 



The source function for pion-produced photons is obtained from 


hv M = 



where (dn/dE p ) is the cosmic-ray proton spectral density. 
We use the approximations 


dE 

dh^ * 2h V s(E„-2hv) , (76) 

dcr / M \ 

dE^^V W \) > (77) 

where cx is the total neutral pion-production cross-section just above 
threshold (~10 -26 cm 2 at 1 GeV). 

These approximations (Equations 76 and 77) and Equation 75 give, 
at hv = 1/2 (M 77 /m)E p , 

% * lftrn o <r .(^)(^) ’ (78) 

where dJ/dE is the spectral intensity of cosmic-ray protons (protons/ 
eV-cm 2 -sec-sr). From Equation 78, note that the number spectrum of 
cosmic pion-generated gamma rays at hv measures the cosmic-ray 
proton spectrum at ~14 (h^) , which is £1 GeV . 
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Recent observations by Clark, Garmire, and Kraushaar (Refer- 
ence 7) with a gamma-ray detector aboard OSO-III indicate that there is 
an appreciable flux of gamma rays £10 2 MeV associated with the ga- 
lactic gas disk. The intensity of these gamma rays is about an order 
of magnitude higher than that which would be expected from the pions 
produced by a cosmic-ray proton beam at several GeV that is similar 
to the one observed near the earth. One possible explanation of the ap- 
parently excessive gamma -ray flux would be to attribute an enhanced 
cosmic-ray intensity to some regions of interstellar space relative to 
what we actually observe near the earth at ~1 to 10 GeV. 

STELLAR STATISTICS FOR X-RAY SOURCES 

In observational stellar statistics, we are concerned with the num- 
ber N* of stars, their apparent intensities I, and the distribution func- 
tion g defined by 

_d / dN* \ 

g = dl \da) ’ (79) 

where fi denotes solid angle. Spatial statistical models again concern 
the number N* of stars; however, here we refer to their intrinsic 
luminosities L and the distribution function p* in space defined by 

, _ A (&A 

P - dL \ d » x ) ’ (80) 


where x is the vector stellar position. The two distribution functions 
of stellar statistics are related by 



g(I, n) dl dn 



(L, R, fl) R 2 dL dR dfi , 


(81) 


where R is the distance from the source to the observer and d 3 x 
= R 2 dRdfi. The parameters L and I are related by the "point source" 
condition 


L 


4ttR 2 I . 


(82) 
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From Equations 81 and 82, one may obtain the relation between g 
and p* as 


g( I > 0) 



L, /L/4tt I , fi)dL 


2(4tt) 3/2 I 2 - 5 


(83) 


For the situation where p* is independent of R (- ih/ 4/7 I ) , the numer- 
ator of Equation 83 is independent of I , and the expression assumes the 
simple form 


g(I, fl) co 1-2-5 _ (84) 

Consider the total number of stars N R * within a field Afi and within 
a distance R. For this case, we have 



(85) 


where n x is the mean spatial number density of stars within the region 
observed. For the case of a finite R and a lower bound to L, there ex- 
ists a minimum observable intensity 1 0 defined by 


AQ 


f 


g dl 


( 86 ) 


This total number of- stars within the examined region N R is to be com- 
pared with the number of stars Nj* for which the observed intensity 
exceeds I : 


N 



g dl . 


(87) 
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The integrated stellar intensity of all photons l AQ observed within 
An is given by 


I 


&a 



gldl . 


(88) 


From Equations 84, 87, and 88, we obtain the relations 


ih) 


1/2 


3 (l 3/2 Nj*) 




(89) 


K) 


1/3 - 


X AQ 


si (n;) 


2/3 


(90) 


From Equations 85 and 90, we get an expression for n, as 


- - 3 (^ g ) 3 

(AO) R 3 (i 3 Nj*) ’ (91) 

In summary, from the observables N * , I, and I AQ , we can obtain 
I fl and N R directly from Equations 89 and 90, respectively. With some 
additional information for estimating R, we can utilize Equation 91 to 
estimate the spatial density of those stars which are defined by the di- 
rection of observation and which substantially contribute to the ob- 
served photon flux within the bandwidth utilized. 

We use the X-ray data of Friedman et al. (Reference 8) to obtain 
I 0 , N r * and n t for two regions of importance to X-ray astronomy, 
Cygnus and Sagittarius. 

For Cygnus, 

Nj* = 8, for I = 0.15 photon/cm 2 -sec, 

I AO ~ 5 photons/cm 2 -sec, for AQ = 0.23 sr. 
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Using Equations 89 and 90 for these Cygnus observations gives 

I Q = 0.08 photon/cm 2 -sec; N R * % 21. 

From the angular disposition of the X-ray sources in Cygnus with re- 
spect to the galactic stellar distribution, Friedman et al. conclude that 
R % 4 x 10 21 cm for this complex of sources. Using this estimate, we 
use Equation 91 to evaluate n x ~ 5 * 10" 63 /cm 3 for the X-ray stars in 
Cygnus. 

For Sagittarius, 

Nj* = 15, for I = 0.3 photon/cm 2 * -sec, 
l ia = 18, for A ,Q % 0.18 sr. 

Friedman et al. estimate that R ^ 8 x 10 21 cm for the Sagittarius 
sources. Using Equation 91, this yields n x ~ 3 x 10~ 63 /cm 3 for the 
number density of X-ray stars in this region. 

The X-ray star densities for Cygnus and Sagittarius appear to be 
quite comparable. In general, this analysis indicates that the number 
density of X-ray stars within the galaxy is <10~ 6 of the number density 
of optical stars. This implies <10 s X-ray stars within the entire 
galaxy. 

This discussion on stellar statistics for X-ray sources is the out- 
growth of some ideas suggested by Professor W. Kraushaar of the Uni- 
versity of Wisconsin. 
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COSMIC PHOTONS 

Appendix A 

USEFUL CONSTANTS 


Electromagnetic Interactions: 


Q * 

r = — - = 3 x 10" 13 cm (Classical Electron Radius) 

e me 2 


= 'J 17 ( r e ) 2 " 7x 10” 25 cm 2 (Thompson Cross-section) 


Quantum Localization of Charge: 


e 2 1 

a = he = 137 (F ine Structure Constant) 


Thermal : 

Boltzmann 


k = 8.6 x 10" 5 eV/°K 

( 1 . 4 x 10" 16 ergs/°K ) 


Stefan 


cr = 5.7 x 10 5 ergs/cm 2 sec (°K) 4 


Photon Energy: 


hv (eV) 


12.4 x 10" 7 /X. (meters) 
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V. PLASMA ASPECTS OF COSMIC RAYS* 

D. Wentz el 

University of Maryland 
College Park , Maryland 

From the point of view of plasma physics, the cosmic rays consti- 
tute an extremely hot and tenuous plasma which coexists in space with 
the "interstellar gas." This concept, which has been stressed by E. N. 
Parker, leads to several interesting consequences. Of particular in- 
terest here is the dynamic role of cosmic rays in determining certain 
properties of the galaxy. We begin with a very simple picture of a 
plane-stratified galactic disk, symmetric about the central plane (z = 0) 
with all physical quantities approaching zero as | z | - <». A gravitational 
attraction toward z = 0 is provided by the stars in the disk, and op- 
posing forces which keep the interstellar gas above the disk (to a scale 
height of order 100 pc = 10 20-5 cm) are provided by its own pressure, 
by the pressure B 2 /8^ of the magnetic field, and by the cosmic rays. 

The field B is thought of for the moment as being parallel to the ga- 
lactic disk and provides a coupling whereby the gas and the cosmic 
rays may interact. The pressures (or energy densities) of the gas, the 
cosmic rays, and the magnetic field are all typically of order 1 eV/cm 3 . 

The interstellar gas is observed to be largely condensed into 
"clouds,” and it is easy to picture an instability that would cause this. 

If one part of a magnetic field line dips slightly in toward the galactic 
plane, the ionized gas tends to flow "downhill" into this region (Fig- 
ure 1). The added weight then further depresses this part of the line, 



Figure 1-Distortion of the mognetic field. 


•Notes taken by D. Hall. 
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while the nearby higher parts now have been unloaded and thus re- 
leased, except for their ends, which are anchored in the gas clouds. 
These "loops" of unloaded field extend well above the galactic plane. 
Now, since the cosmic-ray "gas" has a "sound speed" comparable to 
the speed of light, and gravitational attraction of this "gas" toward 
z = 0 is completely negligible, one would expect it to fill up the mag- 
netic field loops very quickly until a uniform pressure is reached. 

But some of the loops must extend so high that their field strength at 
the top is insufficient to contain this cosmic-ray pressure. These 
loops are thus "inflated" indefinitely, providing escape routes through 
which cosmic rays can leave the galaxy. In this case there would exist 
a gradient of cosmic-ray pressure directly related to the time scale 
of the escape process. At any rate, the galactic disk cannot be in a 
state of hydrostatic equilibrium. 

An attractive feature of this picture is that it provides a regulating 
mechanism to keep cosmic-ray and magnetic pressures comparable in 
the disk: too few cosmic rays will be unable to inflate the field very 
well, and will collect; too many will inflate the field rapidly, and let 
themselves out. 

It is now clearly important to investigate how fast the cosmic rays 
escape. First, observational evidence (the abundances of light nuclei) 
indicates that those we see have been in the disk for an average of 
about 10 6 years. Searches for sidereal diurnal variation have found 
anisotropies $1 percent even at high energies ("TO 6 GeV), indicating a 
net streaming velocity less than 600 km/sec and, therefore, a distance 
from source to escape of only 'TjOO pc. This has been a puzzle: on the 
one hand, the observed field seems to lie along the spiral arm, and 
escape along the spiral arm would mean paths »600 pc; on the other 
hand it is difficult to understand how the particles could diffuse across 
the field so efficiently as to reach the surface of the disk before trav- 
eling 600 pc. However, J. R. Jokipii and E. N. Parker now propose 
that the observations of mean fields parallel to the galactic plane are 
compatible with a model using "stochastic lines of force" (References 1 
and 2). They believe that with reasonable step sizes in the random 
walk a typical line of force anywhere in the galaxy could reach the 
surface somewhere within about 600 pc of its course along the spiral 
arm. 

Now we are faced with a paradox: the magnetic field line on which 
we lie should be connected to a cosmic-ray escape route only about 
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2000 light years away, and cosmic rays should be escaping by blowing 
up the inflated fields at speed c; this fact should be communicated to 
us with a time lag of only 2000. years, so that we should observe the 
cosmic rays streaming toward their exit. But we do not. The reason 
is probably that the streaming is limited to a speed much less than c 
by some sort of dynamical friction, and the remainder of this chapter 
concerns the question of how this friction may arise. 


The interstellar medium constitutes a magnetized plasma which 
will support hydromagnetic waves; the Alfven velocity characterizing 
these waves is v ft ~ 10“ 4 c. Of particular interest are those waves 
that propagate in the direction very nearly along the undisturbed mag- 
netic field, since the others are subject to certain strong damping 
mechanisms. The two transverse modes propagating along the field 
are properly classified by their circular polarization direction. Be- 
cause v A «c, these waves are practically stationary from the point 
of view of a relativistic cosmic-ray particle. Therefore, we may 
picture the wave-particle interaction with the aid of the diagram in 
Figure 2. This stationary circularly polarized wave superposed on a 
static field B 0 has a magnetic field which exerts a force (q/c) vxBj 
on a particle of velocity v; this force has a nonzero time average if 
and only if the gyroperiod for the particle in 
the field B 0 is equal to the time it takes the 
particle to travel one wavelength along the 
field line. That is, there is a resonance I do-' 

condition 


kv 


ii 


n 


q B o 

TTnc 


IN 


for strong interaction; particles of mo- 
mentum p and pitch angle 8 ~ cos 1 u rel- 
ative toB 0 interact with those waves with 
wavenumber 


0 


UT 


B j (wave) 


qB 0 qB 0 
k “ cp n ~ cp/x 

Thus a positively charged resonant par- 
ticle, traveling upward with phase such that 


Figure 2-Diagram of station- 
ary circularly polarized wave, 
the wave field spiraling in 
the sense ofa proton with in- 
dicated v n . 
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the component of its velocity perpendicular to B 0 is into the paper at 
those points of Figure 2 where the arrows (i.e., Bj) point to the right, 
is decelerated, but one with velocity out of the paper at the same 
points is accelerated. If all phases are equally likely, there is, to 
lowest order, no net change in the total momentum of all resonant 
particles taken together. But consider particles nearly in resonance: 
those going slightly too fast with "phase in" have their resonance im- 
proved, and so do those that are too slow and have "phase out." ("Fast 
and out" or "slow and in" particles become less resonant and have neg- 
ligible effect.) If there are more fast particles than slow ones, there 
is a net deceleration of particles; therefore, the wave will gain mo- 
mentum. (Remember that these statements concern acceleration 
parallel to B 0 ; in a static magnetic field, the energy of the particles 
can never change, so the process considered here is one of scattering, 
primarily. Energy changes are slower by a factor v ft /c .) 

With a spectrum of waves present, a particular particle inter- 
acts with different waves from time to time, having v„ sometimes in- 
creased and sometimes decreased. This scattering in pitch angle 
tends to keep the distribution of cosmic-ray momenta isotropic; it can 
largely destroy the anisotropy that would be generated by a gradient 
in cosmic-ray density, and therefore keep the net streaming velocity 
low. 


We now outline a quantitative calculation of the effects described 
above.* 


First, the so-called quasi-linear theory of plasma fluctuations 
provides an equation for the behavior of the expectation value of the 
particle distribution function f(r, p, p., t) in the presence of a spec- 
trum of waves of the type discussed above: 


ii 

dt + v 


Vf 


27 t 2 e 2 d 

M 1 

(■*.) 

1 1-M 2 1 

(dj_ df\ 

c d/J. 

1 

p 2 


\dfl + c P dp J 


♦For more detail see D. Wentzel, Astropbys. ]. 156, April 1969. 
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Here M is defined by requiring that the total energy density in the wave 
magnetic field be 


f 

•'n 


M(k) dk . 


Some terms of higher order in v A /c have been omitted. 

Second, an equation is written for the wave amplitudes: 


1_ <9M(k) 

M(k) dt 


C ! 



J d/xp 2 dp (1 - fj. 2 ) 


8 






i 


'2 r 3 
B 0 


k 2 M(k) 


Here C x and C 2 are calculable constants, n e is the number density of 
the background plasma, and n c r the number density of the cosmic rays, 
with f normalized to unity. The first term describes the reaction upon 
the waves when they scatter the particles; clearly if di/d/x is sufficiently 
large and positive, this gives wave growth. The second term repre- 
sents wave damping by the particular process that we believe is most 
important here: the nonlinear conversion of Alfven waves into sound 
waves, which are then quickly damped. 

Third, we hypothesize a cosmic-ray density gradient (ultimately 
maintained by a source and a sink which are not shown explicitly) that 
can be represented by a scale height which may be a function of energy, 
L( p/mc ) E . We also use observational knowledge to specify that the 
main part of the cosmic rays is described by the isotropic distribution 
function f . = ip" 7 with y - 4.5. Then we ask for a simultaneous 
steady- state solution of these nonlinear equations: 


0 , 


0 . 


<9f 

dt 


<9M 

dt 
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It is found to be 



M = M 0 (k/k 0 ) - ( 7 “'>'~e)/2 . 


The interpretation of the second term within the brackets is that the 
cosmic rays may have streaming velocities up toyv A /3 - 1.5 v ft ~50 
km/ sec without causing instability and building up an appreciable wave 
spectrum; larger streaming velocities must be accompanied (and con- 
trolled) by waves. For permissible values of g, the anisotropy is an 
increasing function of energy. 

This solution seems to confirm some of the ideas we have put 
forward about the physical processes. Cosmic rays of sufficiently 
high energy may flow readily from one place to another and have ap- 
preciable anisotropy. But the pressure of the cosmic rays is due 
mainly to those rays of energy below 10 GeV; and these are prevented 
(by waves they generate themselves) from streaming much faster than 
v A . Hence, the rate of filling and rising of magnetic loops above the 
galactic disk is definitely limited. In fact, as B 0 becomes weaker, the 
instability and scattering become stronger and slow down the excess 
streaming more and more (perhaps until the wave energy becomes 
comparable to B 0 2 /8v, at which point the calculation is invalid). There- 
fore, we may have an even stronger self-regulating mechanism than 
was described earlier. 

It should be possible to investigate certain other interesting ques- 
tions in much the same way. One example is energetic electrons; as 
they synchrotron- radiate, they tend to become anisotropic (although 
with no net streaming). Again, the waves so generated will limit the 
anisotropy. Preliminary calculations indicate that strong anisotropy 
may be found in the local cosmic-ray electrons when it becomes pos- 
sible to extehd our observations to energies of order 10 4 GeV. The 
same question is of some interest in connection with radio galaxies. 
Another problem along these lines is whether, and for how long, par- 
ticles can become trapped in magnetic bottle configurations. 
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VI. RELATIVISTIC STOCHASTIC PROCESSES 

J. R. Wayland 

Department of Physics and Astronomy 
University of Maryland 
College Park, Maryland 


INTRODUCTION 

In high-energy astrophysics there are many phenomena which in- 
volve stochastic processes. One of the most often studied is that of 
cosmic-ray origin. To understand the concepts involved, we will de- 
velop the classical theory of random processes along the lines set forth 
by Wang and Uhlenbeck (Reference 1). We will then follow a parallel 
path to arrive at the relativistically correct formulation of the theory. 

THE CLASSICAL THEORY OF STOCHASTIC PROCESSES 

The process x( t ) is said to be random in time if the variable x 
does not depend in a definite way on the dependent variable time ( t ). 
This means simply that if you observe the process at different times, 
you find that the process is described by different functions. We can 
define the process x( t ) by 


4> x (xt) dx = probability of 
time t; 

^2 ( X 1 t V X 2 t 2 ) dx i dx 2 = 


^3 

If we look at an xt plot 
the distribution at different 


the following set of probability distributions: 


finding x between x and x + dx at 


probability of finding x between x t 
andxj+dxj at time t j , and between (lb) 
x 2 and x 2 + dx 2 at time t 2 ; 

(lc) 

of the trajectories of events, we can sketch 
times as in Figure 1. Because these are 
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probability distributions, it follows that 

<£> > 0 ; ( 2 &) 

4> n (xj tji x 2 t 2 • • • X n t n ) is a function symmetric in the set of 

variables x x , x 2 , • • • , x n but not (2b) 

necessarily in the set x x tj, • • • , x n t n ; 

<£ k ( X 1 x k t U ) = | dx k + l ‘kn^n ( x l t l "■ X n t n) * (2«) 

The condition expressed by Equation 2b results from the require- 
ment that 4> n is a joint probability. Equation 2c arises from the re- 
quirement that each <f> n must imply all previous 4> k when k < n . 

If the cause of fluctuations is independent of time, the process is 
said to be stationary in time. Then in Figure 1, the distribution func- 
tions would aH be the same at the times 1 x , t 2 • • • . The probability 
distributions in Equations 1 now become: 

4> l (x) dx = probability of finding x in the range (x, x + dx) ; 

<£ 2 (x x x 2 t) dx t dx 2 = joint probability of finding a pair of values x 
within the ranges (x x , x x + dx x ) and (x 2 , 
x 2 + dx 2 ) at times that are an interval 
t = 1 2 - t x apart; 


etc. 
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We can assert from the foregoing equations that averaging over an 
ensemble will give us the time average. 

We can use the conditions given in Equations 1 to classify random 
processes. A purely random process is one in which the successive 
values of x are not correlated; i.e. , 

$2 ( X ! V X 2 t a ) = *1 ( X 1 tr) ( x 2 t 2 )’ etC - 

Thus in terms of an xt plot, the distribution function at various values 
of x and t is completely independent of other values of x and t. Then the 
distribution function at each value of x and t will be ( xt ) . This is a 
meani n gful concept as long as t is discrete; however, when t is con- 
tinuous, x x and x 2 will be correlated for t sufficiently small. 

By the conditional probability , P 2 (x 2 | x 2 , t), we mean the proba- 
bility that given x : one finds x in the range (x 2 , x 2 + dx 2 ) after an 
interval of time t. Let us now consider stationary processes only. 

(In the notation used, a bar separates the variables that are given from 
those for which the probability is to be found.) We note the following 
properties: 


^ 2 ( x i x 2 t ) = M x i) P 2 



(Xj | x 2 , t) ; (3a) 

0 ; (3b) 

= l ; (3c) 

| x 2 , t) dxj . (3d) 


A Markov processe s one for which the conditional probability that 
x lies in the interval (x n , x n + dx n ) at time t nJ given that x is equal to 
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x i» x 2 > ’ ’ * > x n-i attimes t v t 2 , ■ • . , t n _ 1 (t n > t nM > . • • t a > t x ), 
depends only on the value of x at the previous time t n _ x ; i.e., 

P (x, t,, ••• x t , | x t ) - P„ lx t .lx t ] . (4) 

In other words, we are considering processes which "remember" 
practically nothing about what has happened in the past. In fact, the 
last event is just on the point of being "forgotten." However, these 
probabilities must still satisfy Equation 3d and also the condition 

LimP 2 (x x I x 2 t) = ^ (x 2 ) . ( 5 ) 


This last requirement simply implies that there are no recurring 

chains. 

Consider two events (x x , t = 0) and (x 2 , t). 


X 



We can write the conditional probability that each will reach some 
intermediate state xt 0 . Then the probability of (x x , t = 0) ->(x, t Q ) 
h>(x 2 , t) is just the product 


P 2 (x x | x, t 0 ) P 2 (x | x 2 , t- t 0 ) . 
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Integrating over all intermediate states x gives 



Xi t o) ^2 ( x I X 2 ’ t *0 ) 


(6) 


This is the Smoluehowski equation. 

Let us now rewrite Equation 6 so that 


P(x 


t + 


At) 


Jdz 2 p (x | Z 2 , t) p(z 2 I y, At) 


(7) 


We are interested in the moments of the spatial coordinate changes in 
a time At, which are defined by 


a n (z, At) 


I 


dy(y - z) n P(z 


y. At) . 


The assumption that one usually makes at this point is that only the 
first and second moments are proportional to At. This allows us, in 
the limit as At - 0, to consider only 

( z > At) 


A(z) 

- Lim At 

At-o 

(8a) 

B(z) 

a 2 (z, At) 
- Lim At 

(8b) 


In most of the problems we will consider, one can actually prove that 
this is the case. We will not give the details here, but from Equations 
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7 and 8 one can derive the Fokker-Planck equation (References land 2): 

<9P a , -i 1 a 2 r -i 

at " “ a y [A(y) p ] + [b(y) p] , ( 9 a) 

for one-dimensional problems; and 

it = ■ [ A i (y) p ] + < 9b > 

1 k,T 

for n-dimensional problems. 


THE RELATIVISTIC THEORY OF STOCHASTIC PROCESSES 

To understand completely relativistic stochastic processes re- 
quires a sound grounding in the fundamentals of special relativity and 
relativistic statistical mechanics. We cannot possibly begin to cover 
the necessary background that this would involve in the short space 
that is available to us. We will instead draw a strong parallel to the 
above classical theory and compare the results of one to the other. 
We quote the results of the relativistic formulation without giving a 
detailed derivation. By taking this approach, we hope that the re- 
sults will be more meaningful in a physical sense. For more com- 
plete discussion of the mathematical detail, see the recent work of 
R. Hakin (Reference 3). 

Let us first define certain necessary quantities. We will con- 
sider a one-particle phase space, if you wish, a Ai- space. To form a 
consistent, but not necessarily unique, set of definitions, we will take 
the metric tensor in Minkowski space-time as 


g°° = "g 11 = +1 

= 0 for m / v 

c = 1 . 


> r 

Mi v 




0, 1, 2, 3 

1, 2, 3 


(10) 


We want to consider a particle whose motion is random in the M-space 


M 


M 4 X U 4 , 


( 11 ) 



RELATIVISTIC STOCHASTIC PROCESSES 


99 


where M 4 is the Minkowski space-time and U 4 is the four- velocity- 
space. We are interested in an arbitrary coordinate system in 
/j-space; thus we will need to consider the system 

(X^, U v ) = X A (A = 1 , ••• , 8) . 

What we are actually concerned with are the trajectories in /i - space 
of a random point, and these are required to be timelike. However, 
in four-velocity space, the trajectories are spacelike. A seven- 
dimensional hypersurface in u-space is spacelike when its normal is 
a timelike eight- vector. Let 2 7 be a spacelike hypersurface imbedded 
in /-a- space such that it cuts all possible timelike trajectories of the 
process in question in u- space. Also let A be a Lebesque measurable 
subset of 2 7 . We want to define a probability 


P(X A e Ad 7 C/a) 


for all A. Then the P's are determined by the distributions <f > i (x A ) 
(analogous to the 4> x (xt) in the classical case): 


P(X A e A C 2 7 C/i) = (x A )d S B , (12) 

J Ac2 7 C> 


where dS B is the differential element of hyperspace, and the V B is an 
eight-vector. The v B are the "tangents" to the trajectories of the 
process. Analogous to the requirement of Equation 2a in the classical 
use, we here require 


4> x (X A ) > 0 , 
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with the normalization 


P(X A e I 7 C/x) 


(, 


'*i ( xA ) d2 s 


= + 1 


(13) 


for all A. 

We can now define the joint probability of finding X * within 
Aj c Sj 7 and X 2 A within A 2 cl 2 7 by 


•[( 


X A £ Aj C 


V)' 


i (X A eA 2 c2 2 7 


)] 


I 


'A 1 xA 2 cS 1 7 xS 2 7 c fl 2 




(X A , X A ) dS B 


(14) 


(a = 1 • • • 8, B. = 1 • • • 16, i = 1, 2) , 


where dX B is the differential element of hypersurface correspond- 
ing to a 14-dimensional subsurface embedded in a 16- dimensional 
manifold, and £ B ! B 2 serves the same purpose as r, B above. This 
corresponds roughly to condition lb. Then from 4> 2 (X A , X A ) we can 
find all the probabilities 

p[(x A , X 2 a ) eAcX 14 Cm 2 ] • 


We have considered the cases of and P 2 separately. Now we 
can expand our definitions to the general case. This will parallel the 
classical definitions given by 1 and 2. 

A relativistic stochastic process in /a- space will be determined 
completely when all the probabilities 

n (x A e A. cX. 7 Cm) . 


p 
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ke positive integer for all 2. 7 are given. Then, as before, we assume 
that all P's are given by the densities 4 k (X A • • • X k A ). Again we must 
require that 

( x i A - > K) > 0 : (15a.) 


4> k (Xj*, ••• X k A ) is not necessarily a symmetric function in (15b) 

Xi A ; 




(15c) 


for all 2 7k and B i = 1 • • • 8k, i =l***k,A = l***8; and 


(V--1C) 


L 


2 7 < n-k >C> (n_k > 


•b„_. 


k ^n(X A ---X A ) dZ B 




(15d) 


for all S 7 ( n-k ) and n > k. Here 


2 7 (n — k) CAt (n-k) 


means a surface, 2" ( n-k ) , is embedded in a space spanned by 
[^nA> * ^(n-k ) a] ' 

The conditional probability is defined by 


^n + l(*0 


y 

B’ a 1B’ 


’ ^b) ’ ^nB I ^0b) ' ' ' ^b) 

(16) 


(compare with Equation 3a). This is the probability density for reach- 
ing the state X 0B by passing through the states X 1B , • • • , X nB . (Note 
that we have reversed the order on the index.) Thus we have a partially 
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ordered sequence of states 


^nB < ' < X 1B X 0B 


In the relativistic case a Markov process must satisfy 


P n ( X 1B’ ’ X nB 1 X 0 b) P 2 ( X 1B I X 0b) 


(17) 


(compare with Equation 4). Again we are led to a Smoluchowski equa- 
tion. In the relativistic case , we find 

P 2 ( X 0 B I X 1b) = ^ A ( X 2 b) P 2 ( X 0bI X 2b) • ( X 2B t X !b) dS A C X 2 B ) ’ 


where X 0B <X 2B <X 1B (compare Equation 6). This leads directly to 
the relativistic Fokker- Planck equation, 


<9 

(9X a 


PA ( X B ) P 2 


i a 

+ 2 <9X b 


° AB ( X b) P 2 


= 0 


(19) 


(compare to Equation 9). Here B A and D AB are the first and second 
moments in the limits of At - 0. 

Under more complete analysis, a number of interesting points 
about relativistic stochastic processes emerge, e.g.: 

1. One can always think of a stochastic process in space as an 
idealization of an underlying dynamic problem if suitable assumptions 
about the scales of the variables are made. 

2. A Markov process may be a consistent idealization if the in- 
variant parameter characterizing the spatial extent of the system, X, 
and the correlation time of the process, T, are related by T » \/c . 
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3. Causality requirements in Minkowski space-time require that 
the Fokker- Planck equation (Equation 19) be applied to a system of 
particles; one cannot apply this form to a single particle. 

4. A relativistic random process must have carefully specified 
tensorial properties, and one should probably also specify how it is 
observed and measured. 

APPLICATION TO COSMIC RADIATION 

Because of the difficulties encountered in trying to obtain an ex- 
act solution to the relativistic Fokker-Planck equation, we will con- 
sider the more often used non-relativistic Equation 9. When con- 
sidering cosmic rays, one normally is concerned with how to obtain 
the very high energies observed. If we consider the case in which 
the variables are time, t, and energy, E, we find that 



where n is the particle density. The first term on the right-hand side 
describes the mean statistical energy change of cosmic-ray particles. 
The second term takes into account the statistical fluctuations in the 
energy change. Normally we can simplify Equation 20 by letting 

= aE , (2la) 


(ae) 

At 


and 


((AE) 2 ) : 


At 


2bE 2 • 


(21b) 
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Thus we must solve 


dn <9n dn 

■37 + an + aE gg ~ 2bn - 4bE gg - bE 2 


d 2 n 
<9E 2 


0 . 


( 22 ) 


The form of the energy operator suggests a power-law solution. There- 
fore, we will apply a Mellin transformation with respect to the energy, 
E. This gives us 


dg 

at 


(s - 1) ag 


(s - 1) (s - 2) bg 


0 , 


(23) 


where 


E s ~ 1 n(E) dE . (24) 


Then we have that 

g = c exp |[(s - 1) a + (s - 1) (s - 2) b] tj . (25) 

The simplest injection one can assume is a delta function in 
energy as 

% s ( E _E o) S { t " *<>) - (26) 

which in the transformed s space becomes 


g(s) - 




^oPo 5 ' 1 S ( t_ t o) • 


(27) 
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The solution for g above must be equal to (27) when t t 0 . Thus we 
find 

g = q Q p 0 s_1 exp |[(s - 1) a + (s - 1) (s - 2) b] rj , (28) 

where r = t - t 0 . The inverse Mellin transformation is given by 

/.c+i^ 

n(E, t) = 2^1 8( s ' OE-'ds . (29) 

J c -i<D 


Dr. F. Jones has pointed out that this integral can be found analytically. 
The result is given by 


n(E, t) 


E 0 i^Trhr 


i*o' y 


pj exp _ (a-b) 2 /4br 


(30) 


where 


- -i-i (*\ + 3b --* . 
y - 4b r Ie o) 2b 

Then we have a power law spectrum in which the exponent is a function 
of time. When br « 1 the spectrum is very steep and is determined by 
the -tn(E/E 0 ) term. If br ~ 1 both terms are important and the result- 
ing spectrum becomes steeper at higher energies. When br » 1 the 
spectrum is determined by the (3b - a)/2b term. If the observed 
spectrum is the result of br >> 1 type sources, then a * -2.2b for 
7 = 2.6 (the observed value). We can only obtain a satisfactory fit to 
the experimental data for the integral spectrum when we assume a 
continuous deceleration (i.e., a < 0). This could be the result of the 
expansion of the source region. Note that we have a balancing of 
continuous acceleration against deceleration. While deceleration 
dominates, an acceleration process must be present. In the case of 
dominating deceleration, the primary cosmic-ray spectrum is en- 
tirely attributable to the consequences of statistical fluctuations. 
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Let us now consider how we can obtain the Fokker-Planck coef- 
ficients a and b. The collision probability is proportional to the 
relative velocity v r of the scattering center (velocity V), and the 
velocity of the particle (velocity v) ; thus 

v r _ (f3 2 + B 2 ~ 2/3B cos 9 - fi 2 B 2 sin 2 d) 1/2 

c 1 ~ ft B cos ’ (31) 


where 8 is the angle between v and V, /? = cv , and B = Vc. it is also 
proportional to the velocity distribution f ( V) of the scattering cen- 
ters, The collision probability </> is then given by 


v r f (V) dVd/x 

'A dVd/x = — ’ (32) 

JJ v r f(V) dVdp. 

where m = cos 8. We will assume that m is isotropically distributed. 

From the foregoing, we will consider the case in which decelera- 
tion can dominate. Let us restrict ourselves, for the sake of illustra- 
tion, to the case of magnetically turbulent scattering centers that are 
receding from each other. Please note that this is not the only pos- 
sible case. The foregoing results are of a very general nature. 

In the case of a spherical expansion from a common center, the 
scattering centers will recede from each other with a velocity kV e /R. 
Here V e is the expansion velocity, k is the mean free path between 
centers, and R is the radius of expansion. The average energy loss 
caused by expansion alone is - (v V e k/c 2 R) E. Let k = kV e /vR. 

If we consider Fermi acceleration, we can show that the fractional 
energy change per collision is given by 


AE 

E 


T 2 (l + 2/3B cos 6 + B 2 ) - 1 . 


(33) 
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We can expand Equations 31 and 33 in powers of B (keeping terms up 
to order B 2 ) and let A - 1 (i.e., high-energy limit). Then using Equa- 
tion 32, we can calculate a and b as defined in Equations 21a and 21b. 
Here we use 


<(£P» 


II 


(AP) n 4 1 dVd/i 


and 

B 2 = |b 2 f(V) dV . 

Thus we find 

At a - 8/3B 2 - 2k + 2/3 kB 2 , 
At b - k 2 - 4k B 2 + 4/3 B 2 . 


(34a) 

(34b) 


We are interested in the case of a - ab, where a < 0 (i.e., a is nega- 
tive). Then using Equation 34, we find 


B 


2 + aB 

e 


l 

R 



J <2-a) + 2(l + 2a) 


1/2 


(35) 


where B e - V e /c. Normally for a nova or supernova, B e ~ 10~ 2 . We 
can make a very rough estimate of the number of scattering centers 
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by 


N 


S C 


number of scattering centers 


4/3.77 R 3 
4/3 77 't 3 



(36) 


We can, with the aid of Equation 31, obtain a plot of B vs RAt as shown 
in Figure 2. 



Figure 2. 


CONCLUSIONS 

We have obtained a relativistic Fokker-Planck Equation. Although 
we have solved a very particular form of this equation, we should point 
out that if one chooses a different set of coordinates, the resulting 
equation will be one of the many forms studied in recent years under 
the heading of relativistic stochastic processes. 
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VII. PROPAGATION OF SOLAR COSMIC RAYS 
IN THE INTERPLANETARY MAGNETIC FIELD 

E. C. Roelof* 

NASA Goddard Space Flight Center 
Greenbelt, Maryland 

The goal of the study of solar cosmic rays is the understanding of 
their production. At present, however, attention is focused on their propa- 
gation in the interplanetary plasma because their production spectrum 
is severely distorted by the time the particles arrive at the earth. The 
interplanetary magnetic field, carried out into space by the expanding 
corona (the "solar wind"), is the main source of this distortion, since 
its irregularities tend to randomize the trajectories of the particles. In 
this messy situation one must be grateful for the fact that the solar 
cosmic rays are so much more energetic than the solar wind particles 
(whose kinetic temperature is less than 10 6 °K, their average energy 
being about 100 eV) that their trajectories may be analyzed as individ- 
ual particle orbits, neglecting cooperative effects in which they would 
modify the fields through which they pass. 

The first step in understanding the propagation of solar cosmic 
rays is to understand the interplanetary magnetic field. Because of 
the huge conductivity of the tenuous interplanetary plasma (<10 
particles/cm 3 at 1 AU), we can speak of the interplanetary magnetic 
field as being "convected" out from the sun because it is "frozen" into 
the plasma. This concept arises from considering the change of flux 
through any closed contour, L, moving with the plasma (at the plasma 
velocity u). This change is due both to (d/d t ) B and to the change in 
the area, since the shape of L will change if u is a function of position. 
The total time rate of change of the flux is thus 
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hence 

d 

at ® = 


Nonrelativistically, the electric field in the moving frame is E' = E 
+ u/cx B. If the plasma obeys Ohm's law, J' - ex' E' . Since con- 
ductivity is_a property of the medium, a' = a; and since the plasma 
is neutral, J' = J + p u = J, resulting in the relation 

J = ct(e + \ u x b) . 


da • (- cVx E) + 0 (-dtxu) • B 
s J L 


: j> dl • (e + | ux b) . 


If_the current is to remain finite while a-a> j we must have E = - u/c 
x B. This is called a "polarization field" because it can be thought of 
as arising in the following way. The Lorentz force per unit charge on 
a particle moving on the average with the plasma is (since (v) = u), 


(f/q) - ( v/c x B } = u/cx B . 


This force tends to separate positive and negative charges and polarize 
the plasma, inducing an electric field E p . The condition that there be 
finite current for infinite conductivity requires E = - u/cx B. (Note 
that a particle not moving with the plasma responds only to its own 
Lorentz force v/cxB, and to the polarization field 2 , the sum of 
which is not zero on the average.) 

The result of this modest excursion into plasma dynamics is that 
d<t>/dt = 0. Thus field lines threading L at t must also thread it at 
another time t 2 . Obviously, great simplification of the "frozen- in" 
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concept is that the magnetic field can be 
deduced from the plasma velocity field. 

The solar plasma corotates with the sun 
while it is close to the sun's surface. 

However, when u becomes comparable 
to Alfven velocity B/V(4v p), the Maxwell stresses cease to dominate 
the mechanical momentum; mechanical angular momentum is then con- 
served, and the plasma escapes radially 
at large distances from the sun. The 
solar wind accelerates rapidly to ~400 
km/sec while the Alfven velocity de- SUN 
creases so that the transition to radial 
flow occurs within ~0.1 AU. 

Problem: Assuming that the plasma flows radially out from the 
sun, where the field B° (0) is B (r Q , ® ), and © is measured from the 
north ecliptic pole, show that the equation of the magnetic field lines 
is the Archimedean spiral <p = - fir sin ©/u and that the field at a 
distance r from the sun can be written 

„ , . . fir sin 0 \ 

B(r, 0) - B°(0) (r 0 /r) 2 (e r n e 0 ) 

(which satisfies V • B = 0) where r Q = the radius of the sun, fi = the 
angular rate of rotation of the sun, and e r and are unit vectors in 
the radial and azimuthal directions, respectively. 

Solution: The Archimedean spiral configuration of the IP field 
line is the locus of the plasma outflow (which is radial, Ar = uAt) 
from a single area element on the rotating source (the plasma moves 
through a central angle A </> = - fiAt sin 0 relative to the area element 
in time At). Note: the field line must lie on the cone © = constant. 
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All points on the locus thus must satisfy 


dr u 

d0 ~ Qsin© 

Integrating, 

fir sin© 


(choosing r = 0 when 0=0). By flux conservation through a moving 
surface wheno--oo, d$/dt = 0; therefore § lies along the spiral, so 
that it must have the form 


B(r, ©) 


o-(r, ©) ^ 


e + r 


d4> 


dr 


spiral 



- a(e r 


fir sin© 
u 



Clearly this contains no explicit 0-dependence, by azimuthal symmetry. 
For the field to be physical, Maxwell's equation V • § = 0 must hold; 
hence, 


1 ± 

r 2 <9r 


r 2 B 


1 3 . 

r sin© 30 Sln © % 


1 d 
r sin© 30 


1 d 

~2di r2a 


a - 


HQ) . 

r 2 


If B(r Q , ©) = B° (0), then 


fl sin© ^6 (®) ■ 


f(0) 


r Q 2 B r ° (0) 
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Since 


V 


fir sin © 


fi - 3 x 10 1 sec 1 , 

r Q ~ 700,000 km , 


and 


we find that 


u ~ 400 km/sec , 


B 0 ° ~ -0.005B r ° ; 


SO 

B°(0) - y/(B r 0 ) 2 + (B°) 2 - B r °(0) ; 


giving, finally, 


B(r, 0) 



fir sin © \ 

u %) ■ 


By construction, V • B = 0. 

These field lines then corotate with the sun. It is now easy to see 
the effect of the polarization field for energetic particles spiraling 
along the lines. A steady electric field produces a drift velocity 

V d = C g X B/B 2 . 

Substituting E p gives 

V d - c(- u/c x B) x g/B 2 - u - g(B • u) = u d ; 

so the field tends to keep the spiraling particles on the corotating 
field lines by giving them a transverse velocity equal to the transverse 
velocity of the lines. 
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For a solar wind velocity of 400 
km/sec, the field line passing through 
the earth (r = 1.5 x 10 8 km) leaves the 
sun at a western longitude 


4 > 

^ o 


nr 

u 


(3x 10~ 6 /sec ) (l , 5 x IQ 8 km) 
4x 10 2 km/sec 


1.13 = 65° • 


The Q in the equations is 2 v divided by the sidereal period. The sun's 
sidereal period is about 24 days , although the rotation period viewed 
from the earth is about 27 days because of the earth's orbital motion. 
The angle that the field makes with the earth-sun line is 


i p - tan 1 — - tan 1 1.13 ~ 50° • 

The general field at the sun (r 0 = 700,000 km) is ~1 gauss; so this 
simple theory predicts that the field at earth is about 5 x 10” 5 gauss. 

All the foregoing theoretical predictions are borne out by direct 
and indirect physical measurements of the average properties of the 
field. However, there are considerable fluctuations about the average 
values. As we shall see later, thesei small fluctuations strongly effect 
charged particle motion. 

Let us now consider the particle observations themselves. Even 
before there was any solid experimentally based theoretical knowledge 
of the interplanetary magnetic field, cosmic-ray physicists had con- 
cluded, from two different phenomena, that energetic particles undergo 
some diffusion-like transport in the interplanetary medium. 
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First, the time dependence of solar cosmic rays following large 
flares sometimes looked like a diffusion curve. The classical three- 
dimensional diffusion equation 



DV 2 N 


has the solution 


N(t) 


( 477 Dt ) 3/2 


-r 2 /< 


4D t 


Noting that In (t 3/2 i) - (r 2 /4D) (1/t) + constant, one multiplied 

the counting rate I(t) by t 3/2 and plotted this product versus 1/t 
(where t = o at the time of the flare). The result was often a good 
fit to a straight line. (Sometimes t a I, where 1 < a <3, gives a better 
fit than a = 3/2.) Interpreting the slope as -r 2 /4D with r = 1 AU, 
the value of D was found to be on the order of 10 2 2 cm 2 /sec. This 
number seemed reasonable, since its use in the classical relation 
between the mean free path and the diffusion coefficient, D = \v/3, 
yields A.-^lo 12 cm = 0.06 AU. Even though the actual scattering 
process was not known, it was comforting to deduce that there were 
"many" scatterings between the sun and the earth so that the diffusion 
picture was self-consistent. However, it was clear that the process 
had to involve magnetic fields, since the mean free path for Coulomb 
interactions is (7rr B 2 ohr n) -1 ~10 16 /n centimeters which is many AU 
for reasonable values of the interplanetary proton density n . 

The second phenomenon that supported the diffusion model was 
the modulation of galactic cosmic rays. Neutron- monitor counting 
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rates showed a long-term variation with an 11-year period which was 
just opposite that of indicators of solar activity (i sol = sunspot number 
or corona green- line intensity, for example). This variation was ex- 
plained on the basis of the convected magnetic irregularities. The 
galactic cosmic rays must "swim upstream" against the irregularities 
which tend to diffuse them. The net flux, then, is that due to the dif- 
fusion process (-DVn) plus that due to simple convection at the ve- 
locity of the frozen-in field (nu) : 


J - DVn + nu . 


In a quasi- steady state, V • J ~ o. If we invoke spherical symmetry 
in the solar system so that D, n, and u are functions of r only, and also 
assume the direction of u to be radial, then V x J = 0. Therefore J is 
a constant; since J = Je r , the only constant value consistent with 
spherical symmetry is zero. This gives the equation for n: 


Vn/n = u/D , 


which may be integrated along any path (L) from r to r Q to give 







n(?) 

- n Q exp 

-1 dt • u/D 

- n o exp 

r r o 

- 1 dr u/D 


_h 


_ Jr 


The modulating region must terminate at some distance r Q at which 
n = n , the galactic intensity. In the approximation u/D ~ constant, we 
can estimate the distance to the boundary if we assume D does not 
change much. From solar maximum to minimum, the average value 
of the solar wind decreases by about 100 km/sec while the neutron- 
monitor counting rates increase by 20 percent. From the equation 


^ Ar , > ^ Ar , 

-1 = 1 _ exp u' \ u min _ u max ) = irK 
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we obtain Ar ~ +(+0.2) (1022)/10 7 ~ 2 x 10 14 cm ~ 10 AU, a number 
consistent with other estimates of the size of the solar cavity. The 
reader can find detailed discussions of all arguments up to this point, 
as well as many other theoretical concepts relating to cosmic rays, 
in Reference 1. 

To make further progress towards the primary goal of under- 
standing the nature of the solar cosmic rays at the sun, a microscopic 
description of particle motion in disordered magnetic fields is neces- 
sary. The calculations involve considerable arguments and the tech- 
niques of the theory of stochastic processes. However, a great deal 
of this may be circumvented by realizing that the velocity process in 
a spatially disordered magnetic field must be a random walk of the 
velocity vector on a sphere in velocity space. The contribution of the 
electric field - (l/c)u x B to scattering is negligible compared to the 
Lorentz force (1/c) vx B for v » u. (However, the electric' field is 
important for long-term (~days) energy loss in the medium.) This 
justifies assuming that the magnitude of the velocity is constant. 

The problem of a random walk on the surface of a sphere was 
solved beautifully by F. Perrin in 1928. He found the differential 
equation for the time-evolution of the distribution function in spherical 
polar coordinates (a, <p, where ,u = cos 8) to be 


H _ 1 \±_t. „\ Bf , 1 d 2 f 

<9t T |j?M ^ 1 -fjfi aq) 2 


The operator on the right is the "transverse" part of the Laplacian in 
velocity space (i.e., that part not involving v). Its eigenfunctions are 
the well-known spherical harmonics, products of Legendre polynomials 
in /a and cissoids in <p. There is, reasonably, one parameter in the 
equation, the "relaxation time," T. Its role is easily seen from the 
Green's function solution 


CO -v 0 

>. 0.. .> - £ ^ .-«*•** £ M**™. 


Z=c 
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lim f 0 (Ai, 9. t) = 8(/i"M 0 ) 8(<P-<P 0 ) • 


For times t » T, an arbitrary initial velocity distribution on the sphere 
will have relaxed to near isotropy. 

For magnetic scattering, a Fokker- Planck analysis leads to the 
Perrin equation if the mean field B Q vanishes. If B Q does not vanish, 
its presence introduces two new ter ms, .yielding the equation 


dt 

dt 


n. 


ii + 1 

<9q> T 


An explicit expression 
particle orbits: 


dfJ, 


dp ( 1_ ^ 2 ) 


<9 f 


1 d 2 f 


1 -p* <9<p 2 


+ 4 


d 2 f 

<9cp 2 . 


for T follows from the statistical analysis of 


T 


2. /me \ 2 v 

17 U ' P i ( k o) ' 





while 


P ll (0)- p i( k n) 
P x( k o) ’ 


where P x (k) andp^ (k) are one-dimensional (measured along the 
mean field) spatial power spectra of a random field component per- 
pendicular and parallel to the mean field, respectively. Clearly the 
new terms vanish for B Q = 0. This equation is valid only in the limit 
of small-angle scattering, which is equivalent to 


M/v « l , 
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where 


n 



and l is the "correlation length" of the random field. The terms "mean" 
and "random" are properly defined relative to the scale of the root- 
mean-square gyroradius, v/0; thus l for low energy particles may be 
much shorter-than l for those of high energy, since the former smoothly 
follow irregularities of scale length much larger than their gyro- 
radius, while the latter see these same irregularities as random per- 
turbations. Thus Z really depends on energy for a given field 
configuration. 

The physical consequences of the additional terms (which vanish 
when B fl = 0) are (1) a simple tendency to gyrate in the mean field with 
frequency , and (2) an enhanced azimuthal diffusion on the velocity 
sphere due to the mean field. Both terms can inhibit spatial diffusion 
across field lines, as we shall see. 

Description of the process in velocity space is not enough— we 
need the phase space distribution function W(r, v, t). Since the velocity 
scattering is small-angle, the change in position, Ar, over a time in- 
terval At « T (the appropriate time scale of the complete equation) is 
~ vAt. Thus the full equation is obtained by replacing d/dt with 


If the mean field is constant (we shall see later why this restriction 
is important), the full equation is 



where the components of v are ( v sin 6 cos cp, v sin 8 sin <p, v cos 9). 

This equation is quite complicated, but fortunately we need only 
coarse I nf ormation to discuss the general aspects of propagation. 
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Therefore degrade the equation to a diffusion process by assuming that 
W is nearly isotropic in velocity space; i.e., we take only the zero and 
first-order terms of a spherical harmonic expansion: 


w = ^ + 4^r (yi -M 2 cos 9 + /r 


sincp J 2 J 


,)• 


where N and J are functions of r, t , and only the magnitude of the ve- 
locity. Substitution in th e equa tions, multipli cation through by the 
orthogonal functions 1 , i\ - /j. 2 coscp, VT~ sincp, and /j., and integra- 
tion over - 1 < fi £. 1 and 0 £ cp < 277 gives the four equations (written 
in vector form): 


dN 

dt 


+ V 


d v 2 

at J + T 


VN 


2 

T 


—~x Cl x ( Cl x T ^ - 

n 2 x o ' o J / 


n 0 x 


Consistent with the quasi-static diffusion approximation, we assume 



0 


so the vector equation may be solved for J. 

j - - (d„ v (| n + d, v I N ) + xn o xVN , 


where 


D„ = v 2 T/ 6 , 


D, 


2(2 + 0 

(2 + O 2 +(fi 0 T) 




2 “II 


B o B o • V ' 


V, 


V - v„ . 
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The last term is like a J x B term and (for B () constant) contributes 
nothing to the diffusion process when J is substituted into the equation 
of continuity 


at = D iVN + D,j V ,, 2 N . 

If B 0 were not constant, the cross term would have introduced the well 
known curvature and gradient drifts . The constraining of the particles 
to field lines is quantitatively shown by the ratio Dj/Dy , which varies 
as (n Q T)" 2 for T » | 2 + £ | . As the relaxation time grows much larger 
than a gyro-period, the transverse scattering becomes negligible. 

The parallel diffusion coefficient can now be expressed in terms 
of the transverse power spectrum of the field 



_v_ R 2 
3v P x (B 0 /R ) ’ 


R - mvc/ q . 


This form suggests the introduction of a mean free path, X, via the 
traditional kinetic theory relation D = Xv/3 . We have also introduced 
the particle "rigidity,” R, which determines the gyroradius, R/B. Note 
that X = R 2 /P i (b q /r) is a function of particle rigidity, while Dy de- 
pends on velocity as well as rigidity. The observed power spectrum 
near 1 AU has roughly the form of a power law, k” m , with m between 1 
and 2. This predicts X « R( 2_m ) . Solar flare proton events above 20 
MeV have been fit with rigidity- independent X, while modulation of 
galactic cosmic rays appears to require X«R. Thus, to within one 
power of R, the theory is consistent with observations . The measured 
values of P(k) also give the right order of magnitude for Dy. 

More information than just flux time histories has been obtained 
from experiments; in particular, the directional flux distribution has 
been measured (usually in the ecliptic plane). To interpret these ob- 
servations, we must return to the full equation and, even more important, 
must allow B fl to be a function of position (e.g., the Archimedian spiral 
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field we previously discussed). In this case the axes in velocity space 
depend on position, since /x is measured from B Q ; therefore, if we wish 
an equation for W(r , v, /x, cp, t), we must take account of the implicit 

dependence of /x and cp on r. Thus the 
total particle derivative with respect 
to r with v held constant involves the 
explicit determination of Vw as well as 
the implicit terms: 



d 

df w 


VW + 


a \ aw 

B o) a/x 


v(v • e ij'j 1 - (v • § 0 ) 2 )(il^) 


aw 

a cp 


since /x - v 


B 0 and v 


ej = Yl - fi 2 cos cp, if cp is measured from some 


vector e x normal to B Q . ^The complete equation becomes, after using 


the identity V(v • S Q ) = (v • V)B q + vx ( Vx fi Q ), 


aw 


aw 1 


dT + v • VW - fl 0 ^+T 


( v • (v • V)B 


v aw cot 
o) a^ j _ 


A. 2 \ dW + 1 

dM ) djj- + \y- 

(v(v • V)e ± ) 


+ 4 


a 2 w 
a cp 2 . 


cot cp aw 
2 acp 


sin cp Y l ~ fJ. 2 


aw 

acp 


Rather than deal with this entire equation, we can simplify it for solar 
flare events, using the following argument. Applying the relation 
given earlier for T in terms of D n and v, we can estimate that T~ 6 D/v 2 
~/3 -2 (minutes) where fi = v/c. At 1 AU, 0. Q ~ 1/2 sec -1 , so T~30/3 -2 
thus (fl Q T) 2 is a large number, so D^D,^ « 1 and nonrelativistic par- 
ticles must follow field lines. This suggests that W must be independent 
of cp , for lines connecting to the flare region, and that the only spatial 
parameter is the distance out along the field line, x. Setting <9W/dcp = 0 
and averaging the complete equation over cp, we have an equation for 
the new function w(x, v, /x, t): 


dvi 3W 

It + ^ dZ 


v , „ v dw 1 <9 , „ \ dw 

dl l 1 ~ /x ) Jjl + T d/J. i l V- ) dfJ- 
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where L - - B d /dB o y&x. The new term arises from the manipulation 


j.r 

277 

* n 


• (V • v)l o = \ (i -M a ) v • l 0 + f^ 2 t 0 


(8 0 -v)B 0 


V / 

2L ( 1_ 


A 2 ). 


since (b q • v)b q is perpendicular to B 0 and V • B q - - (l/B 0 ) <9B 0/ /dx. 

The physical significance of the new term becomes clear if we 
allow T -* <», in which case scattering becomes unimportant. Then we 
are left with a linear first-order partial differential equation which is 
immediately soluble in terms of the constants of integration of the 
characteristic equation 


dx 

p-v 



We can find /x(x) immediately by cross- multiplying and substituting 


l 

L 


JL^o 

Bq 


2m4u + (1 -M 2 ) ^ dB 0 = 0 , 

which integrates to (l - /j. 2 )/b q , constant. This is the well-known 
first adiabatic invariant (the magnetic moment M), so called because 
the product of the current and the area of a charge in helical motion 
is 


M - q(statcoul) • (qB /277 me) ( sec" • 77(11^ c/qB ) 2 (cm 2 ) 

= (c/q) (mv 2 /2) (1-/7. 2 ) B . 
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Thus this term represents the tendency for the M of a particle to re- 
main constant along a field line (dl^/dx = 0), so that in a diverging 
field, (like the interplanetary field), the particle tends to "run out" 
along the field lines. Only the presence of velocity scattering (T <®>) 

can prevent this strong collimation. For 
example, if the pitch angle 6 = tt/2 for 
a particle at the sun (B 0 ~ 1 gauss), at 
1 AU (B 0 ~ 5 x 10' 5 gauss) the pitch 
angle would be 

0 - sin0 = /Bearth/Bsun = °-°° 7 = °‘ 4 ° 

This competition between collimation 
and scattering is regulated by the single 
parameter a - vT/2L, as may be seen by multiplying the transport 
equation through by T 



d w (?w 

r Tt + ^ Js 


- a(l -p. 2 ) 


dw 

djJ. 
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where ds = dx/vT. When a-0, scattering clearly predominates and 
we have random transport; but when a -*°o there is deterministic mo- 
tion (M = constant). 

Problem: Show that e' 2as e a ^ is a solution and sketch the function. 
Solution: Substitution leaves 

da 

O' 2s) gi - 0 , 


so if a = constant (or if s da/d s ~ 0), this function is a solution to 
the steady state problem. This is indeed the dominant mode (s > 1) 
for steady outflow from the sun. The simplest assumption about 
propagation would be s ~x/vT (T - constant); then as ~ V2L ~ 1, 
independent of particle properties. Note, however, that the above 
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assumption is consistent with the condition on the solution (s da/d s 
= 0) only for x » vT, since s d a/d s = x 3 a/d x ~ vT/x. 

To study time-dependent outflow, it is again simplest to degrade 
the equation with a weak anisotropy approximation. 

Problem: Let 


n(x 


t) 3 . 

2 + 2v K x ’ 


t > 


and obtain the continuity and current equations by multiplying the w 
equation by 1 and v/i and integrating - 1 <fi < 1. In the quasi- stationary 
approximation, dj/dt = 0, obtain the diffusion equation. What are 
these equations for B = B Q ( r 0 /r) 2 e r ? 

Solution: 


dn <9 j _ 1 . 

Jt + Jx ~ ~ L J ’ 


d j v 2 dn 2_ , 

at + 3 9x T J • 


Since - 1/L = (1/B) dB/dx, the first equation can be written 


§T + B i/B = 0 . 


In a curvilinear coordinate system defined by the field lines, such that 
one coordinate (x) is the distance along a line, it can be shown that for 
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a vector j which is always parallel to B, 

v ■ T = B — j/B - 


This result depends on V • B vanishing. Combining the density and 
current equation in the diffusion approximation, we have 

dn _ d 1 (v 2 T\dn 
d t B <9 xB\6 ) dx ’ 

where the (position-dependent) diffusion coefficient is identified in 
the operator V • (DVn), For a radial field we recover the familiar 
equations of isotropic diffusion 


in 1 a 

<?t r 2 d r r 


dj vf dn_ 
dt + 3 ar 


<9n _ 1 d 

d t r 2 d r 
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dn 

dr 
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Two aspects of propagation have not been considered here. First, 
electric fields have been neglected. This is justifiable for energetic 
solar cosmic rays since the characteristic time for appreciable energy 
loss is on the order of the duration of events (or longer) ; however, the 
effect on galactic cosmic rays, which can dwell in the interplane- 
tary field for days, is of considerable importance. Second, the one- 
dimensional equations appropriate to the description of solar events 



(i) 
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do not include any information on the relative positions of the field 
lines. The feet of the field lines undergo some disordered convection 
in the photosphere due to the turbulent velocity fields there. Thus 
field lines originating at positions 1 and 2 may later be at positions 1 
and 2'. If the solar wind has convected these fields out continually, it 
is clear that they will be quite tangled in interplanetary space. Another, 
and probably more important, mechanism leading to disordering of the 
spiral pattern operates in the Alfven region, where the outgoing coronal 
plasma makes its transition from corotating to radial outflow. In fact, 
the transport of the field from the photosphere through the chromosphere 
and the corona to the Alfven region is far from understood. 

Now that we have covered the basic ideas in particle propagation, 
let us turn once more to the observations. The problem of galactic 
modulation is still unresolved and a discussion of its present status 
would be beyond the scope of this discussion. Suffice it to say that 
the diffusion-convection theory gives a generally acceptable explana- 
tion of modulation. 

With respect to solar cosmic rays, we need only up-date the re- 
view by C. E. Fichtel and F. B. McDonald (Reference 2). (Note: The 
a +/3 in Equation 18, p. 391 of their paper should simply be a.) As 
solar (and experimental) activity increased after 1966, the last period 
covered in the review article, events became more frequent and more 
complex. The role of low-energy particles (protons <20 MeV, elec- 
trons >40 keV) has increased in importance, ironically, because they 
do not usually propagate in the manner of classical diffusion. They 
scatter considerably less, as is evidenced, for instance by the prompt 
arrival of 40 keV electrons about 30 minutes after the optical flare. 

This corresponds to almost direct spiraling out the field line, since 
yS> 0.4 for these electrons, and line-of-sight transit is 8 minutes for 
light. Low energy particles also exhibit strong pitch-angle anisotropies, 
again showing that there is much less scattering than at higher en- 
ergies. There is also a strong statistical tendency for 40 keV elec- 
trons to be present after flares on the western half of the sun's disc 
(i.e., longitudes between 30 and 90 degrees W), but not after flares on 
the eastern half. Thus low-energy particles serve as tracers ' of the 
interplanetary field lines. Moreover, they are often long-duration 
tracers, since they can be continually produced by the sun for a num- 
ber of days. These ideas arose from the analysis of recurrent low- 
energy events in which enhanced fluxes appear on successive solar 
rotations about 27 days apart. One well documented nonrecurrent 
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event gave a composite sample of all the foregoing modes of propaga- 
tion (Reference 5). By means of three satellites at different positions: 
IMP III, and OGO III near the earth, and Explorer XXXIII near the moon 
(Figure 1), the following picture of this complex event was deduced: 

The high energy protons (Figure 2a-II) arrived in a clean, classical 
diffusion mode after the flare at 0027 UT, July 7, 1966. However, 15 
MeV protons (Figure 2a-IH) show a hint of a second maximum a day 
after the flare. This tendency is clearer in the protons sensed by 
OGO HI (Figure 2b), in which diffusion dominates at 32 MeV while the 
second-day maximum is practically all that is seen at 3 MeV. The 
low-energy electrons (Figure 2a-IV and V) show prompt arrival, but 
also a later maximum with an even sharper peak superimposed (about 
0800 to 1200, July 8). The interpretation (Figure 3) is that there was 
prompt diffusive propagation of all components at early times, but 


TO SUN 



Figure 1. Locations of the IMP-111, Explorer XXX 111, and OGO III 
satellites during July 7 to 9, 1966 projected onto the ecliptic 
plane. OGO III is at high geomagnetic latitude (14 to 37 degrees) 
and therefore outside the magnetosphere during most of this 
period. 
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Figure 2. Part (a)shows: ( 1 ), the major solar flares and radio bursts, geomagnetic 
disturbances and K p for the period July 7 to 9, 1966; (II), (III), (IV), and (V), the 
counting rates of particle detectors on IMP-1 1 1 and Explorer XXXIII, Part (b)shows 
the count rates of six different energy channels of the 0G0- 1 1 1 detector for the 
same time interval. 
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Figure 3. A schematic of the core and halo of the energetic particle 
fluxes on July 7. The core and halo are convected by the solar wind 
past the earth. The lower left graph indicates the spatial variation of 
the particle fluxes. The upper diagram shows the thin diffusing layer 
around the sun postulated by Reid in 1964 (Reference 3) and Axford in 
1965 (Reference 4). p is the distance away from the flare measured in 
the diffusing layer. 


that the sun continued to produce low-energy particles through July 8 . 
Since the flare was at longitude 48 degrees W, the time for the flare 
region field lines to corotate to a position connecting to earth was 
about (60° - 48°)/(13° per day) ~ l day. These field lines were filled 
with electrons (>40 keV) and protons (~10 MeV) so that the detectors 
actually recorded the spatial profile of the particle fluxes. 
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Problem : Using the contents of this paper, explain the propaga- 
tion characteristics of the four groups of solar particles discussed. 

Solution : For protons above 20 MeV we know there is diffusion 
with D ~ 10 22 cm 2 /sec. Since /3 > 0.2 and the scale length is r/n for a 
field falling off as r~ n , the anisotropy parameter is 

2L c /3L ~ ) (0 2 ) (i s xio 13 ) 0-3 ’ 


so there is only weak velocity anisotropy, consistent with classical 
diffusion. Despite the fact that constraint to field lines is considerable 
since 
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the particles can still have a time history compatible with classical 
three-dimensional diffusion, because for B® r" 2 the one -dimensional 
diffusion equation is identical to the equation for isotropic classical 
diffusion 
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dt 


~ D 
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Furthermore, if (k) cc k“ 2 , then X. is independent of rigidity, as in 
classical transport theory. Thus energetic protons appear to diffuse 
classically. 


Electrons of energy 3 to 12 MeV have a much smaller gyroradius 


than energetic protons: for example R 
However if, again, P t ( k ) a k“ 2 , then 
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and MeV electrons should diffuse like MeV protons. Since their values 
of fi are of the same order, the electrons also would not have large 
pitch angle anisotropies. 

On the other hand, as we go to lower-energy protons and electrons, 
less scattering and higher anisotropies are observed. Since the smaller 
velocities tend to decrease D and a, the only way D and a can increase 
is for Tto increase. A strong anisotropy of a = 1 (actually, a~ 5 has 
been observed for protons and electrons) gives vT - 2L ~ 1 AU so that 
the distance traveled in a relaxation time is comparable to our distance 
from the source. Clearly the diffusion approximation cannot be valid, 
although the one-dimensional equation (with steady-state solution 
e - 2 as e a^ p ro bably is still a reasonable description of this low-energy 
propagation. 

We must end this discussion of the propagation of solar cosmic 
rays with a strong qualification. Another explanation of at least the 
early portions of high energy flare particle events is a very real pos- 
sibility, although it is not incompatible with the ideas we have already 
presented. It is quite possible that flare particles spend a good part 
of their time after acceleration in the outer solar corona (Reference 2) 
before streaming out into interplanetary space (Reference 3). The 
time histories predicted by this "two-region" model are almost indis- 
tinguishable from those predicted by pure interplanetary diffusion. 

Only better measurements than we have now of velocity anisotropies 
throughout events can resolve the question because effective interplane- 
tary diffusion implies small anisotropies (a < 1) while the predominance 
of' coronal storage with weak interplanetary scattering implies large 
anisotropies (a > 1) . 
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VIII. INTERPLANETARY DUST 
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R. A. Schmidt has defined interplanetary dust as "microscopic par- 
ticles occurring outside the limits of the earth's atmosphere." The 
particles are variously called cosmic dust, cosmic spherules, microme- 
teorites, micrometeoroids, nanometeorites, meteoric dust, meteor itic 
dust, primordial dust, zodiacal dust, interstellar dust and galactic dust 
(References 1 and 2). In Schmidt's definition, all reference to the phys- 
ical properties and origin of the dust is avoided. An interplanetary dust 
particle is too small to be detected as a meteor when it enters the 
earth's atmosphere: the radius is less than 1 mm. But it can be de- 
tected by its impact on a satellite or space probe and, together with 
many other particles, it causes the zodiacal light. However, the "typ- 
ical" zodiacal particle may be different from the particle typically de- 
tected by its impact on a micrometeoroid satellite (Reference 3) . 

Interplanetary dust has been a subject of interest for several 
centuries (Reference 4), and particularly in this age of space flights 
and interplanetary experiments— because dust is a hazard for such 
enterprises: the life time of delicate satellite instruments is limited 
because of erosion caused by interplanetary dust. To estimate the 
erosion rate we need to know the size distribution, mass per unit vol- 
ume, and velocity distribution of the particles. These quantities have 
not yet been determined very accurately, but an upper limit to the space 
hazard has been obtained (References 5 and 6). A second reason for 
our interest in interplanetary dust is that it is a kind of "blanket" which 
blocks the diffuse star light from our view. Of particular interest in 
this connection is the total mass of dust in a column vertical to the 
ecliptic plane. Interplanetary dust is of further interest to us because 
it may be an important source of certain elements in the earth's crust. 
Interplanetary dust particles are an integral part of our solar system, 
and our knowledge of their physical properties, origin and evolution 
significantly increases our understanding of the solar system . 
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EXPERIMENTAL EVIDENCE ABOUT INTERPLANETARY DUST 
The Zodiacal Light 

There is a faint cone of light centered near the ecliptic plane and 
extending from the solar limb to the antisolar point (References 7 
through 10). This is the zodiacal light which is sunlight scattered by 
dust particles between the sun and the earth and in regions beyond the 
earth's orbit. Near the anti solar point the zodiacal light brightness is 
slightly enhanced: This effect, called gegenschein has been variously 
attributed to (Reference 11): 

1. An increase in particle concentration along the extended sun- 
earth line (a dust tail of the earth) 

2. A localized dust cloud at a libration point on the extended sun- 
earth line (References 12 and 13) 

3. An enhancement in the differential scattering cross section of 
dust particles near 180 degrees (glory) (References 14 through 18) 

4. Continuum emission from a gas tail of the earth (analogous to 
a comet tail) (References 19 through 21) excited by the solar wind. 

The third explanation is satisfactory, for the following reasons: 

(1) a dust tail of the earth requires a large geocentric dust cloud which 
we do not believe exists. (2) The dust concentration at a libration point 
required to explain the brightness of the gegenschein is dynamically im- 
possible (References 22 and 23). Also, the center of the gegenschein 
should always lie exactly on the extended sun-earth line; this is not the 
case (Reference 24). The emission from a gas tail should be enhanced 
following a solar flare, but this is not observed (Reference 22). 

It has been suggested that the zodiacal light— all or part of it— is 
caused by a geocentric dust cloud (References 25 and 26), but this is 
unlikely to be the case because: 

1. The suggestion is based on the rates of particle detection by 
satellite-born acoustic sensors (References 27 through 29), and these 
rates are believed to be spurious (Reference 30). 
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2. No process has been discovered by which such a cloud is 
created (References 31 and 32). 

3. The manner in which the zodiacal light brightness varies with 
ecliptic latitude is in conflict with the hypothesis (Reference 33). 

The zodiacal light has been studied not only from the ground but 
also with rockets (Reference 34), balloons (References 8 and 35 through 
37), and from orbiting spacecraft (References 36 through 39). In prin- 
ciple, it is possible to determine the average physical characteristics 
and the size spectrum and spatial distribution of zodiacal dust from the 
observed brightness (as function of elongation and ecliptic latitude), 
polarization, color, and spectral features of the zodiacal light. A single 
dust model which explains all the observations satisfactorily has not 
yet been found. The relatively large polarization of the zodiacal light 
(up to 25 percent) (Reference 40) presents a special problem: Although 
dust models have been constructed (References 41 through 44) which 
explain the polarization (in these models the average particle radius is 
less than 1 micron and the size distribution, given by a power law— 
dN(s) = Cs' p ds, has a relatively large spectral index p> 4) such models 
are in conflict with results based on satellite impact measurements 
(Reference 45) . The satellite data agree much better with a model in 
which the average particle radius is several microns, andp < 3. The 
polarization may also be caused by the peculiar physical shape of the 
particles. We know yet very little about the shape and the surface de- 
tails of zodiacal dust; until recently, the theory of the zodiacal light 
had been based on the assumption that the particles are spherical. Now, 
calculations of scattering functions are made for non-spherical par- 
ticles (References 44 through 49). Estimates of the bond albedo of 
zodiacal particles range from 0.005 to 2:0.6 (References 7, 10, 16, and 
50). 


The concentration of zodiacal dust in the ecliptic plane at 1 AU is 
between one and several hundred particles per km 3 , and the mass/volume 
between -10 -24 gm/cm 3 and -10~ 21 gm/cm 3 . It is assumed that the 
concentration varies with distance r from the sun as r"- 5 , and fi is 
probably less than 3* (References 7, 10, 34, 35, 41, 42, 44, 45, and 51). 
The fact that the zodiacal light cone appears to be fairly narrow sug- 
gests that the dust particles are strongly concentrated toward a plane 


Powell et al. (Reference 47) have recently put forward the idea of distinct dust belts 
such as a Venus-Earth belt, a Mars-Earth belt, etc. 
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close to the ecliptic. However, the zodiacal light brightness depends 
on both the concentration and the scattering function of the particles; it 
is a quantity integrated along the line of sight and the strong concen- 
tration therefore, may be an illustion. After a study of the isophotes, 
some investigators (References 45 and 51) concluded that the average 
orbital inclination of the particle orbits may be as much as 30 degrees. 
This is contrary to the assumption of a thin dust layer. On the other 
hand, recent measurements (Reference 52) of the Doppler shift of 
Fraunhofer lines in the zodiacal light can be understood better in terms 
of small orbital inclinations (Reference 53) . We do not know whether 
the dust is distributed symmetrically about the ecliptic plane, the in- 
variable plane (off by 1.5 degrees), or the plane containing the solar 
rotational equator (off by 7 degrees). Evidence for each possibility has 
been cited (References 10, 54, 55, and 56). Symmetry properties of 
the dust distribution are directly related to the origin of dust or to the 
effects of perturbations on dust orbits. 


Satellite Impact Experiments 

Special sensors have detected the impact of interplanetary dust on 
satellites and space probes. The sensors are usually one of two types: 
microphone detectors (acoustic detectors, sounding boards) and pene- 
tration detectors (pressure cells, grids, wire cards, parallel plate 
condensors). In the first type, the particle is detected by its sound of 
impact; in the second, by its puncturing a thin film, grid, or plate. 
Neither detector gives any information other than that the particle has 
sufficient mass or velocity to overcome the detector threshold. De- 
tectors are also being developed which measure the velocity or mass 
of an impacting particle (References 57 through 61). 

The accurate calibration of acoustic and penetration detectors has 
proved difficult. In the United States, acoustic detectors are considered 
to be momentum-sensitive but are considered energy sensitive in the 
U.S.S.R. (Reference 60). Penetration detectors are energy sensitive, 
but there is no agreement concerning the proper functional relation be- 
tween particle mass and velocity and detector threshold (References 62 
and 63). In a micrometeroid impact experiment one does not measure 
the dust concentration , but measures instead a kind of flux. An inter- 
pretation of the counting rate in terms of the distribution of particle 
radius and velocity requires some assumption about these distributions. 
These assumptions are usually based on studies of the zodiacal light. 
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Although our only knowledge about the impacting particles is that 
they overcome the detector threshold, the satellite experiments give us 
quite a bit of information about interplanetary dust. This is because we 
can determine how the counting rate varies with (1) satellite altitude 
(distance from the earth), (2) detector threshold, and (3) type of de- 
tector. The early acoustic rates obtained near the earth were much 
higher than expected on the basis of our present knowledge about zodi- 
acal dust (References 27 through 29), and they were interpreted in 
terms of a very high dust concentration near the earth* (Reference 28). 
Attempts to explain why such a dust belt should exist have failed (Ref- 
erences 25, 31, and 64 through 71). It was also discovered that the un- 
expectedly high counting rate near earth could be related to the thres- 
hold requirement of the detectors (References 45, 65, and 72). But even 
so, a very steep size distribution (and an unreasonably low geocentric 
velocity— which is the speed of dust relative to earth at 1 AU) would be 
required (References 31 and 45). Also, the variation of the acoustic 
rates with threshold was compatible only with a flat size spectrum 
(p < 3) (References 45 and 73). Eventually, the acoustic rates were 
identified with thermal noise in the detectors (References 30 and 75), 
and this put the dust belt to rest. 

The penetration rates near the earth, and the acoustic and penetra- 
tion rates far from the earth, agree fairly well with rates predicted for 
models of the zodiacal dust such as given in Reference 7 by van de Hulst 
(p = 2.6; s > lfj.\ mass/volume = 3 x 10“ 21 gm/em 3 at 1 AU) and for a 
geocentric velocity of 5 to 15 km/sec (References 45, 73, 76, and 77). 

In some cases the impact rate on the earth's morning side was dif- 
ferent from that on the evening side, and there were seasonal varia- 
tions as well (References 78 through 81). From such asymmetries one 
can derive valuable information about the sizes of interplanetary dust 
particles and their orbits (References 45, 73, and 82). Consider, for 
instance, dust particles which are partly supported by solar radiation 
pressure: their orbital velocity is smaller than the Keplerian velocity. 
Hence if they are in nearly the same orbit as the earth, the earth over- 
takes them and the influx of this dust into the atmosphere is greater on 
the morning side than on the evening side. The satellite impact rate, 
however, is not necessarily greater on the morning side since particle 
trajectories merge behind the earth. A systematic analysis of the sa- 
tellite data with respect to time variations has not yet been undertaken, 
but the necessary theory has been developed (Reference 73). 

‘Although there was no convincing evidence from the altitude dependence of the impact 
rates (Reference 74). 
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Dust in the Atmosphere and in Sediments 

Since the end of the last century, dust has been variously collected 
with balloons, airplanes, and rockets from the atmosphere and from 
ice and snow cores (References 3, 4, 83, 84 and 85). In the atmosphere, 
dust is collected on specially prepared surfaces which are exposed for 
a predetermined interval of time and altitude range. Prior to sending 
them aloft, the surfaces are carefully cleaned of laboratory dust and 
other contaminants. Some surfaces are purposely not exposed to the 
atmosphere and some others are kept in the laboratory; these serve 
as controls. Such experiments require the utmost cleanliness in the 
laboratory in order to avoid contamination. Nevertheless, in some ex- 
periments the control panels were found actually to contain more par- 
ticles than the panels exposed in the atmosphere! (Reference 85). We 
do not yet possess reliable criteria by which to identify a particle as 
either extraterrestrial or terrestrial. But in recent very carefully 
conducted rocket collection experiments (Reference 87), nearly all par- 
ticles were much like terrestrial dust. Some suggestions about the 
possible composition of interplanetary dust are based on the composi- 
tion of meteorites which often contain anomalous amounts of nickel, 
cobalt, manganese, and sodium (References 85, 88, and 89). From the 
amount of dust collected and classified as extraterrestrial, the rate of 
accretion of such dust by the earth can be deduced. Estimates vary be- 
tween 10 and 10 7 tons/day (Reference 1). 

Interplanetary dust possibly contains considerable abundances of 
certain nuclides depleted in the earth's crust. Analyses of deep sea 
sediments indicated Al 26 in relatively large quantities (References 86, 
90, 91 and 92). In interplanetary dust, Al 26 is produced by solar cosmic 
rays (10 MeV protons) from Mg 26 , Al 27 , and Si 28 . By assuming equi- 
librium saturation of the dust particles with Al 26 (this requires a suf- 
ficiently long exposure to the cosmic ray flux) and by making reasonable 
estimates of the relevant cross sections and of the cosmic ray flux, an 
Al 26 equivalent dust accretion of some 10 3 tons/day by the earth was 
derived (Reference 93). Lately, however, it was suggested in Refer- 
ence 94 that the Al 26 may be entirely of terrestrial origin (from at- 
mospheric argon). Furthermore, no anomalous amounts of Al 26 were 
found in Greenland ice cores (References 95 and 96). The minute con- 
centration of the nuclide in sediments and ice makes precise meas- 
urements of course difficult. The Al 26 theory requires a stony dust 
composition; therefore, only a lower limit to the accretion rate can be 
obtained. Iridium and osmium are also depleted in the earth's crust. 
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From their abundance in deep sea sediments, a dust accretion of less 
than 10 3 tons/day (Reference 97) has been derived. From the abun- 
dance of Cl 36 , an accretion of slightly more than 10 3 tons/day has been 
derived (Reference 98). Other estimates of the dust accretion rate are 
based on the abundance of He 3 , Ar 36 , Ne 21 , and noble metals. (Refer- 
ences 85 and 99). 

In moderate and high latitudes one often detects dust layers, called 
noctilucent clouds, in the upper atmosphere; they sometimes cover an 
area of several 10 6 km 2 , and they are very thin (~2 km) and have a 
well-defined altitude (~80 km). The concentration of dust in these layers 
is between 10” 2 and 1 per cm 3 , and the particle radius is between 0.05 
and 0.5, (References 100 through 102). Noctilucent clouds are also ob- 
served in the antarctica (Reference 103). The dust concentration can- 
not be explained in terms of straightforward sedimentation of dust (of 
interplanetary origin) in the upper atmosphere (Reference 104). Sea- 
sonal variations of the NLC and their peculiar structure, together with 
the fact that the particles are often coated with ice, indicate that the 
clouds are meteorologically conditioned (References 105 through 107). 

Attempts are now being made to determine the dust content of the 
upper atmosphere directly by measuring the return signal of a laser 
beam. Dust layers have been found at altitudes from 60 to 90 km (Ref- 
erences 108 through 110). The dust concentration is of the same order 
as that in noctilucent clouds. The extent of the laser back scattering 
is not known. 


Evidence from Comets, Asteroids, 

Meteors and Meteorites 

Some ideas about the properties of interplanetary dust can also be 
based on our knowledge of dust comets, asteroids, meteors and mete- 
orites. Type II (dust) comets continually lose dust, and at a rate prob- 
ably greatest near perihelion. Dust particles are also created in col- 
lisions between asteroids, and Poynting-Robertson drag (see below, 
DYNAMICS OF INTERPLANETARY DUST) moves them toward the 
Sun. Meteorites and meteors can be taken as evidence for the exis- 
tence of much finer matter in the solar system. The size distribution 
of meteoric particles is steeper than the size distribution in some 
zodiacal dust models. Opik (Reference 112) has suggested that both 
types of particles may have a common origin, but a non-uniformity in 
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the size distribution may arise because large particles are more easily 
disturbed by Jupiter than small particles. But such a non-uniformity 
may not actually exist (Reference 113). It has also been suggested that 
zodiacal dust has a different composition than meteoric particles (Ref- 
erence 111), but there is no direct evidence on this point. Concerning 
composition and shape of interplanetary dust, little is really known. 

(1) The Al 26 theory is based on an assumption that Interplanetary dust 
is stony, and says nothing about other types of dust nor about the shape 
or structure of the stony dust. Some of the features of the zodiacal 
light have been successfully explained (2) with models of purely metallic 
dust and with models of purely dielectric dust, also with mixtures of 
dielectric and metallic particles. (3) From the satellite impact rates 
we cannot yet deduce much about the shape, density, and composition 
of the particles. (4) The composition, structure, and shape of particles 
collected in the atmosphere, deep sea sediments, and snow and ice 
cores are well known; but the extraterrestrial origin of these particles 
is unproven. 


THE DYNAMICS OF INTERPLANETARY DUST 

Several non-gravitational forces affect the orbits of interplanetary 
dust, and the evolution of an orbit depends on their relative importance. 
Some forces are steady over many orbital periods while others change 
rapidly. The physical properties of the dust particles also change with 
time. 


Radiation Pressure and Poynting-Robertson Drag 

A net force, P, is exerted on a dust particle by the solar electro- 
magnetic radiation if the incident photon momentum is not all scattered 
in the forward direction. If E is the energy flux in erg/cm 2 -sec, and c 
the speed of light, then 


P - vs 2 Q pr E/c , (1) 

where vs 2 Q pr is the radiation pressure cross section of the particle, 
and Q pr depends on particle size, composition and shape; it can be larger 
than unity. For various types of spherical particles, Q r has been cal- 
culated from Mie theory (References 114 through 118). Since E is 
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proportional to r -2 , P does not secularly change any orbital parameter; 
however, the orbital velocity is by a factor (1 - P/G) 1/2 smaller than 
the Keplerian velocity where G is the gravitational force on the particle. 
P/G is proportional to (1/s), and very small particles cannot be in el- 
liptic orbits. The value of s fbr which P = G, called radiation pressure 
limit , is approximately 0.2/a. A particle is in a parabolic orbit and it 
would leave the solar system, barring other effects, if its total energy 
is zero. This requires 


P/G ~ 1 " r/ (2a 0 ) , 


( 2 ) 


where a 0 would be the semi-major axis of the orbit if P = 0, (Refer- 
ence 119). The particle radius corresponding to (2) is roughly twice 
the radiation pressure limit. 

The dust particle also experiences a retarding force (drag). Con- 
sequently, its total energy decreases with time and it tends to spiral 
toward the sun (References 120 and 121)— this is the Poynting-Robertson 
effect. (The orbit also precesses, but there are no other changes.) The 
semi-major axis a and the eccentricity e both decrease and are related 
by (Reference 122) 

a(l-e 2 ) e~ 4/5 - constant . (3) 


For a circular orbit (e = 0) the time T in which the orbital radius de- 
creases from & 1 to a 2 is 


T = 710sS a t 2 (l-a 2 2 / ai 2 )Q p V (4) 

years, where s is the particle's radius in p, 8 its density in gm/cm 3 , 
and a 2 and a 2 are in AU. For an initially nearly parabolic orbit with 
aphelion Q[AU] and perihelion q[AU] , the total life time is 


T 


2.3 x 10 3 s 8 QV* q 3/2 q-1 


(5) 
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years. The Poynting-Robertson drift velocity, da/dt, is proportional 
to s' 1 ; the particle mass is proportional to s 3 ; therefore, the net flow 
of dust toward the sun across a unit surface at r is proportional to the 
total surface area of dust in a unit volume at r. One can compute the 
rate at which interplanetary dust falls with the sun directly from the ob- 
served surface brightness of the zodiacal light, and therefore, assump- 
tions about the particle size distribution are not necessary. The rate 
is approximately lton/sec (Reference 123). There is some uncer- 
tainty about this number, as our knowledge of the dust distribution 
above the ecliptic plane is incomplete. 


Planetary Perturbations 

The effect of distant perturbations on dust orbits by the planets is 
insignificant because the life time of a particle is short compared with 
the time scale associated with such perturbations. However, during a 
close encounter with a planet, a large change in the particle's orbital 
velocity can occur. This change takes place instantly as compared 
with an orbital period. In the planetary frame of reference, the total 
energy of the particle is conserved; therefore, the planetocentric speed 
U of the particle before and after the encounter is the same. A second 
encounter with the planet, if it occurs, has the same value of U, and the 
effect of successive encounters is that orientation of the U -vector 
changes randomly. The particle's orbital semi-major axis, eccen- 
tricity and inclination are related by (Reference 111) 

3 r/a - 2^(a/r) (l - / cos i = const. 


(Tisserand's criterion). The probability of an encounter at a planeto- 
centric speed U and impact parameter d is, per orbital revolution (Ref- 
erence 111), 


P 


(d/r) 2 


U ( 
4v v f sin i 


(7) 


where v r is the radial component of the heliocentric velocity of the 
particle. After N encounters, the net change in an orbital parameter is 
proportional to N 1/2 . Since Poynting-Robertson drag removes the 
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particle relatively quickly from the region of space where an encounter 
with the planet is possible, N is relatively small; and therefore, the net 
change in any orbital parameter is small. Under certain conditions, 
however, a particle can be ejected from the solar system by a single 
encounter. 


Corpuscular Radiation Pressure and Coulomb Drag 

Interplanetary dust carries a surface charge which is a conse- 
quence of the absorption of solar UV photons and the accretion of ions 
and electrons from the solar wind. The charge is usually discussed in 
terms of an equilibrium potential V Q , and V Q is a few volts, certainly 
not more than 10, perhaps even slightly negative. V Q is nearly inde- 
pendent of particle radius and varies little with distance from the sun 
(References 124 through 128). 

Solar wind ions and electrons which do not directly hit the dust 
particle nevertheless impart momentum to the dust particle because 
of their Coulomb interaction with the surface change— if the impact 
parameter is less than the shielding distance A. from the dust particle. 
Approximately, \ = 0.5N e ~ 1/2 , where N e is the electron concentration 
in the solar wind. We account for physical impacts and Coulomb inter- 
action by defining the total cross section a. For the i th constituent of 
the solar wind (charge q 4 ) the cross section is given as follows (Ref- 
erence 104): Let V* = V Q let K* be the kinetic energy, T* the mean 
thermal energy of the constituent, u its mean thermal velocity, and w 
the solar wind velocity. Then 


cr/rrs 2 = 1 +V*/K* + (v*/K* ) 2 + ln[(K*/V*) (2>/s)} if K* » T* 

= (2/3) (u/w) [l +V*/f +2 (v*/l*) 2 In{{l*/V*) (2>/s) - 0. 58)j 


if K* « 1* (7) 

The solar wind velocity is 300 to 500 km/ sec (References 129 through 
131), V Q a few volts, and the solar wind temperature 10 4 to 10 6 °K (Ref- 
erence 132). We finder ~ ns 2 . The total force on the dust particle is 
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therefore almost entirely due to impact of ions. It is much smaller 
than the solar radiation pressure. However, the corpuscular drag is 
comparable with Poynting-Robertson drag: 


solar wind drag _ g 

Poynthing-Robertson drag — 1-3 x 10 F , (8) 

where F is the solar wind flux at 1 AU, 10 8 to 10 9 ions/cm 2 -sec (Ref- 
erence 104). The solar wind co-rotates with the sun to some extent 
(References 129 and 130), and therefore the drug on a dust particle is 
somewhat smaller than given by (8). It has been suggested that there 
may be a large number of very small dielectric dust particles in inter- 
planetary space which are transparent enough that they do not contrib- 
ute much to the zodiacal light (They are essentially "invisible."): these 
may be the noctilucent cloud particles. Such particles would experience 
very little Poynting-Robertson drag but they would always experience 
corpuscular drag and would therefore spiral into the sun much like 
under Poynting-Robertson drag— barring other effects, of course. If 
the solar wind were to co-rotate with the sun fully, then the corpuscular 
drag on the zodiacal particles would be strong enough to result in a 
pancake-shape zodiacal cloud (Reference 113). This is not observed. 


Interaction with the Interplanetary Magnetic Field 

The Lorentz force exerted on a dust particle due to the interaction 
of its surface charge with the interplanetary magnetic field is 

L = V Q s(v +w) x B/c , (g) 

where v is the particle's orbital velocity and B the field strength. The 
term involving w arises because the magnetic field is "forzen" into the 
solar wind and is carried along with it (Reference 132). Because of the 
sector structure of the field (Reference 133), the direction of L changes 
by 180 degrees many times during a single orbit of the particle. The 
change in an orbital parameter can therefore be described by a one- 
dimensional random walk (References 125 and 128). Of course, changes 
in the orbital parameters are inter-related; however, if we consider a 
large number of particles with initially similar orbital elements then 
we can calculate the r.m.s. change in any orbital element independently. 
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As an example, consider the inclination after a time T of initially 
circular orbits and having negligible inclinations: An r.m.s. value of 
90 degrees can be expected in (Reference 104) 


1Q10 §2 a 4 a 

t V Q 2 B 2 w 2 (10) 


years, where 8 is in gm/cm 3 , s in p, r in AU, V Q in volts, B in gammas, 
w in km/ sec, and t is the time (in years) which the particle spends in 
one sector of the magnetic field. For typical values S = 3.5, s = 1, 
r = 1, t = 1/90 (4 days), V Q = 5, B = 5 (Reference 131), and w = 500, 
we find T = 7 x 10 4 years. This is much longer than the Poynting- 
Robertson drift time da/dt given by (4). But we notice that T is pro- 
portional to s 4 , whereas da/dt « s 2 : clearly, the Lorentz force is very 
important for sub-micron size dust (s < 1/x). 

Sputtering, Evaporation, Collisions 

The typical energy of solar wind ions is 1 keV (Reference 130); 
therefore they sputter atoms from the dust particle surface, and con- 
sequently the particle is gradually eroded. A recent estimate of the 
erosion rate is ds/dt = 0.05 r -2 A/year (Reference 134), r in AU. 

The time scale associated with sputtering is much longer than the 
Poynting-Robertson time , and there are no significant dynamical ef- 
fects associated with sputtering. 

The temperature of an interplanetary dust particle, T(r), is not 
simply the black-body temperature. For dielectric particles (Si0 2 ), 
Over (Reference 135) finds 


T(r) = 120 r" 1/2 °K; r in AU- (11a) 

Whereas for metallic particles, Becklin and Westphal (Reference 136) 
find 


T(r) 


418 r"°' 27o K; 


r in AU. 


(lib) 
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According to the kinetic theory of gases, the rate of surface mass loss 
by evaporation is proportional to P v T~ 1/2 where P v is the vapor pres- 
sure at a temperature T, and 


l°g 10 P v (T) _ 13.66 - 22200/T for dielectrics (Reference 135) 


12.78 - 20040/T for metals (Reference 136). ^2) 

(For vapor pressure laws for Fe and Ni, see also Opik- (Reference 137). 
A dust particle at r shrinks at the rate 


ds/dt 



(13) 


where m is its molecular weight. A metallic particle at r = 1/50 AU 
(= 4.4 solar radii), ds/dt ~ 1/j. per orbital period. Within a distance of 
a few solar radii from the sun, there are no solid particles. Some 
peculiar dynamical effects are associated with thermal evaporation of 
dust (Reference 128): consider a particle with spirals toward the sun 
under Poynting-Robertson drag; because it evaporates rapidly near the 
sun, its drift speed da/dt increases. Also P/G increases; eventually 
the particle has too much velocity to remain in a circular orbit, then 
the semi-major axis and eccentricity increase. Eventually, the total 
energy becomes positive and the particle leaves the solar system. At 
the distance where the particles stop spiralling toward the sun and 
move out again, we expect an increase in the local dust concentration. 
Curiously enough, some enhanced infrared emission has been detected 
at about 5 solar radii and it is definitely thermal radiation by inter- 
planetary dust (Reference 138). 

Collisions between dust particles affect the size distribution of 
interplanetary dust and the distribution of their orbital elements. As 
a consequence, both distributions vary with distance from the sun. 
Sufficiently small collision fragments are expelled from the solar sys- 
tem by radiation pressure (cf. Eq. (2)). The total amount of dust lost 
from the zodiacal dust cloud in this manner is estimated to be 10 to 
100 tons/sec (References 104, 123, and 139)— i.e. considerably more 
than evaporates near the sun. Collisions are particularly important 
for large zodiacal particles, and the chance that a sub-micron particle 
ever makes a collision is very small (Reference 104). 
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Because of collisions, and also because of the impact of solar wind 
ions, interplanetary dust particles rotate rapidly. We calculate the rate 
of rotation by assuming that the rotational vector performs a random 
walk due to ion impacts. Then the rotation rate after a time t is (Ref- 
erence 104) 


v(t) = 10~ 21 (Ft) 1/2 w (Ss 3 ) sec -1 (14) 

(cgs units). A dynamical effect associated with the rotation of dust is 
the Yarkovsky-Radzievskii effect (Reference 111): A phase lag in the 
surface temperature exists between the morning and evening side of the 
particle and therefore the vapor pressure on these sides is not the 
same. The particle experiences a net force and is accelerated or de- 
celerated in its orbit, depending on the sense of rotation. 


THE ORIGIN OF INTERPLANETARY DUST 

There are several possible sources of interplanetary dust: aster- 
oidal collisions, fragmentation of comets, capture of interstellar dust, 
condensation of interplanetary gas, ejection of lunar surface material 
by meteoroid impact, and a primordial dust cloud residing in the outer 
regions of the solar system. Some of these may not be important. We 
particularly want to know which contribute most to maintaining the 
interplanetary dust cloud observable to us, and what the size and spatial 
distribution of dust from those sources would be. We also want to know 
if the present interplanetary dust cloud is a permanent feature of the 
solar system. 


Asteroidal Collisions 

In order to estimate the rate of production of zodiacal dust in aster- 
oidal collisions, we need to know the rate of collisions as a function of 
asteroidal size, and the total mass and size distribution of debris re- 
sulting from a single collision. Although the collision probabilities for 
given asteroidal orbits are well known (References 110 and 111), the 
collision rate is not . This is because we do not know the distribution 
of asteroids with radii less than 1 km. There is no agreement as to 
whether the asteroid belt is an adequate dust source (References 35, 

104, and 119). Concerning the size distribution of asteroidal dust it 
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has been found that if the distribution can be described by a power law 
N(s) « s'p, then p = 3.8 (Reference 141). The debris drifts toward the 
sun under Poynting-Robertson drag; since the drift velocity is propor- 
tional to ( 1/s ) , the effective size distribution of asteroidal dust has a 
spectral index p = 2.8 (in good agreement with the "flat" dust models 
deduced from zodiacal light observations) (References 7, 16, 46, 113, 
and 142). The velocity of ejecta is much smaller than the orbital ve- 
locity of the asteroids, and therefore the inclination of the dust orbits 
will be asteroidal (mean inclination: 9.5 degrees (Reference 143)). 

The spatial distribution of interplanetary dust of asteroidal origin, as 
determined by Poynting-Robertson drag, is such that the concentration 
is proportional to r - 1 (References 104 and 139) . 

Dust from Comets 

The relation of the rate of emission of dust by comets to the dis- 
tance from the sun is not well known. For comet Arend-Roland an 
emission rate of 75 tons/sec at perihelion is estimated, but it has been 
pointed out that because of solar radiation pressure most dust particles 
of cometary origin do not remain in the solar system (Reference 119). 
The distribution of eccentricity and semi-major axis of cometary dust 
is therefore quite different from the distribution for the comets; how- 
ever, the inclination of the orbits of cometary dust is much like the in- 
clination of cometary orbits (mean: 15 degrees, for short period 
comets; long period comets have a random distribution of inclinations 
(Reference 143) . If a single comet is responsible for the interplanetary 
dust cloud, then (because of Poynting-Robertson drag) the dust concen- 
tration should vary as (a) r' 1 for r <q, where q is the comet's perihelion 
distance from the sun, and (b) r -2 - 5 for r > q (References 104 and 139). 
Probably many comets contribute to the zodiacal cloud, but the con- 
tributions may not be equally important. 


Capture of Interstellar Dust 

The mass/volume of interstellar dust near the solar system is 
- 2 x 10" 26 gm/cm 3 , and the solar motion relative to the nearby 
stars is about v = 20 km/sec (References 143 and 146). (This is also 
the solar motion with respect to interstellar Ca II clouds.) If we take 
20 km/sec to be the solar motion relative to interstellar dust clouds, 
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then the dust mass/volume at a distance r from the sun is 



/ 4.5\ 1/2 

- 2 x 10“ 26 x ( 1 + — p- 1 gm/cm 3 , r in AU , (15) 

where M is the mass of the sun, R its radius. Clearly, a mere gravita- 
tional accumulation of dust near the sun does not account for the ob- 
served zodiacal dust cloud. It has been pointed out (Reference 145) that 
dust can be captured into an elliptic orbit via a single encounter with a 
planet (mainly Jupiter). The capture is permanent because Poynting- 
Robertson drag removes the particles from that region of space where 
a second encounter can occur which could result in the ejection of the 
particle from the solar system. (Even if ejection did occur, a particle 
would have already spent a very long time in the solar system.) The 
typical interstellar dust particle is much smaller than the typical zo- 
diacal dust particle. For the two most seriously discussed interstellar 
dust models, the average radius is 0.16 and 0.3 (References 146 and 
147)., and for these particles P/G— 2. They can never be captured into 
elliptic orbits. But even if we neglect radiation pressure, we find the 
capture rate to be insufficient to explain the zodiacal cloud. An addi- 
tional difficulty concerns the distribution of orbital inclinations of 
captured dust; the distribution is random. However, orbits which have 
a relatively small aphelion are more likely to have a small inclination 
(Reference 148). Such particles contribute most to the dust concentra- 
tion in the inner solar system. Since they constitute only a very small 
fraction of the total amount of captured interstellar dust, it must be 
concluded that interstellar dust is not a major source of supply for the 
zodiacal cloud. 


Lunar Ejecta; Condensation in situ; Primordial Remnants 

The velocity at which debris is ejected from the site of impact of 
a meteoroid on the lunar surface is much smaller than the escape ve- 
locity from the moon's surface (2.2 km/sec); hence little debris es- 
capes from the moon. Furthermore, the earth's orbital velocity is 30 
km/sec, and therefore those particles which do escape are in earth-like 
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orbits and will eventually be accreted by the earth (and by the moon to 
a lesser degree). As a consequence, very little dust of lunar origin 
should be found beyond 1 AU, certainly much less than is observed. 

Formation of solid particles by condensation of the interplanetary 
gas is probably not an important source of dust: The solar wind prevents 
condensation in the inner solar system (sputtering). Farther away from 
the sun the gas density is too low and the condensation time too long. 
Some interesting questions can be asked about the dynamics of con- 
densation nuclei far from the sun: Because of their small size, radia- 
tion pressure and radiation drag are unimportant; corpuscular drag 
may be important but the Lorentz force on these particles is very 
large (cf. Eq. (10)), thus the particles may wander around randomly until 
they encounter a comet or other large body or enter the inner solar 
system where they are quickly eroded by sputtering. 

The sun may be surrounded at large distances by a dust cloud much 
like the (hypothetical) comet cloud. There is no evidence for or against 
this hypothesis. If we consider only Poynting-Robertson drag, then a 
particle now at 1 AU was at a distance of 1.3 x 10 3 /(Ss) 1/2 AU some 
4.5 x 10 9 years ago (cf. Eq. (4)). 
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